The n-body problem in General Relativity up to the second post-Newtonian order 

from perturbative field theory 
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Motivated by experimental probes of general relativity, we adopt methods from perturbative (quan- 
tum) field theory to compute, up to certain integrals, the effective lagrangian for its n-body problem. 
Perturbation theory is performed about a background Minkowski spacetime to 0[(v/c) 4 ] beyond 
Newtonian gravity, where v is the typical speed of these n particles in their center of energy frame. 
For the specific case of the 2 body problem, the major efforts underway to measure gravitational 
waves produced by in-spiraling compact astrophysical binaries require their gravitational interac- 
tions to be computed beyond the currently known 0[(v/c) 7 ]. We argue that such higher order 
post-Newtonian calculations must be automated for these field theoretic methods to be applied 
successfully to achieve this goal. In view of this, we outline an algorithm that would in principle 
generate the relevant Feynman diagrams to an arbitrary order in v/c and take steps to develop the 
necessary software. The Feynman diagrams contributing to the n-body effective action at 0[(v/c) 6 ] 
beyond Newton are derived. 



I. INTRODUCTION AND MOTIVATION 

In this paper we are concerned with the problem 
of describing the gravitational dynamics of arbitrary 
?i > 2 compact non-rotating bodies moving in a 
background Minkowski spacetime. By assuming non- 
relativistic motion, this problem can be approached in 
a perturbative manner, by approximating these com- 
pact objects as point masses and calculating the effec- 
tive lagrangian L e s[{x a , v a , v a , . . . }] for their coordi- 
nates {x a |a = 1,2, ...,n.} and their time derivatives 
{Vai v a , . . . }, up to some given order in the typical 
speed v of these n objects: 1 Newtonian gravity starts 
at O[v ] and the Einstein- Infcld-Hoffman lagrangian 
[l|, that describes the precession of the perihelion of 
elliptical orbits, is of 0[v 2 ] (1 PN). 2 The n body prob- 
lem at 0[v 4 } was first tackled by Ohta et al. Q. Some 
computational and coordinate issues encountered there 
were clarified by Damour and Schafer In the lat- 
ter, some integrals could not be evaluated. A portion 
of these were later performed by Schafer [1], so that 
currently, up to the n = 3 case is known. But to 
know the effective lagrangian for arbitrary n at this 
order, one needs to further calculate the integrals for 
the n = 4 case. As we will see later, once L e ff is known 
up to ?i = 4, the arbitrary n-body lagrangian will fol- 
low from a limited form of superposition. 

We will examine this problem using perturbative 
field theory techniques introduced in [l3j. The mo- 
tivations are two-fold, both of them stemming from 
experimental probes of gravitational physics: one re- 



quires the 2-body effective lagrangian to higher than 
C[tj 7 ], and the other may need the n-body counterpart 
at 0[v% 

Gravitational Waves Detection The recent years 
have seen an array of gravitational wave detectors such 
as GEO, LIGO, TAMA, and VIRGO coming online. 
These experiments seek to detect gravitational waves 
produced by binary black holes and/or neutron stars as 
they spiral towards each other. Within their frequency 
bandwidth, these detectors are able to track the fre- 
quency evolution of the gravitational waves from these 
binaries over O[10 4 ] orbital cycles and hence make very 
accurate measurements. To be able to do so, however, 
theoretical templates need to be constructed so that 
the raw data can be integrated against them to deter- 
mine if there is a significant correlation. Via a gen- 
eralized Kepler's third law relating orbital frequency 
to the binary separation distance, these templates are 
based on energy balance: the rate of energy loss of 
these binaries is equal to the power in the gravita- 
tional radiation emitted. Both the notion of energy 
and expressions for the flux of gravitational radiation 
require the knowledge of the dynamics of these bina- 
ries, which in turn is encapsulated in their effective 
lagrangian. Due to the high accuracy to be attained, 
this effective lagrangian needs to be computed up to 3 
PN and higher. 3 

Currently, the dynamics of compact astrophysical 
binaries is known up to 3.5 PN. 4 (See §1.3 of Blanchet 
Q and the references therein.) To obtain the dynamics 



1 We use units where all speeds or velocities are measured in 
multiples of the speed of light, i.e. c = 1. 

2 The nomenclature is: 0[v 2 Q] Q PN. 



3 Blanchet [f| offers a review of the post-Newtonian framework 
and its relation to gravitational wave experimental observ- 
ables. 

4 The half integer PN order lagrangians, scaling as odd pow- 
ers of v relative to Newtonian gravity, describe dissipation - 
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at 4 PN and beyond is a challenging task. Because of 
the need to regularize the divergences that arise from 
approximating compact objects as point particles, one 
may wish to engage field theoretic methods to han- 
dle them. Such a pursuit was initiated in [l3| . where 
it was shown how to carry out the field theory effec- 
tive lagrangian calculation in a systematic manner by 
first doing some dimensional analysis. One of the main 
thrusts of this present work is to attempt to make as 
methodical as possible such a route in post-Newtonian 
calculations. In particular, we advocate using the com- 
puter to automate the process, so that at the end only 
those Feynman diagrams that truly require human in- 
tervention are left for manual evaluation. Given the 
computational effort required at 4 PN and beyond, we 
believe this is necessary not only to save time and en- 
ergy, but also to reduce human errors. For example, at 
such a high PN order, even the derivation of the neces- 
sary diagrams will itself be non-trivial - the reader not 
convinced of this fact is encouraged to look at appendix 
[Cl containing the 3 PN diagrams - but an efficient im- 
plementation of the algorithm that we will sketch in 
the main body of this work will allow automatic gen- 
eration of Feynman diagrams to arbitrary PN order, 
modulo computing power. 

Solar System Gravity Closer to Earth, the 
Einstein-Infold- Hoffman lagrangian, at 0[v 2 ] beyond 
Newton, is routinely used to compute the solar sys- 
tem ephemerides, and to analyze spacecraft trajec- 
tories and space based gravitational experiments. A 
range of experiments, such as the new lunar ranging 
observatory APOLLO, proposals to land laser rang- 
ing missions on Mars and/or Mercury, and spacecraft 
laboratories - GTDM, LATOR, BEACON, etc. - will 
begin to probe the non-Euclidean nature of the solar 
system's spacctime geometry beyond I PN by measur- 
ing the timing and deflection of light propagation more 
precisely than before. (See, for instance, Turyshev [1] 
for a recent review.) 

Within the point particle approximation, both the 
solar system dynamics and its geometry can be gotten 
at simultaneously by computing from general relativ- 
ity the effective n-body lagrangian. Because general 
relativity is a non-linear field theory, knowledge of the 
2 body lagrangian is not sufficient to deduce its n-body 
counterpart, as superposition is not obeyed. That the 
n-body L e ff encodes not only dynamics {x a [t]} but also 
the geometry g^ [£, x) can be seen by adding a test par- 
ticle to the ro-body system. 5 Denoting the latter's mass 



and coordinate vector as M e and y M = (t, yY respec- 
tively, in the limit as M e tends to zero relative to the 
rest of the other masses in the system, we know its 
exact action has to be 6 



M e / dtxg^ 



dyV dy u 



dt dt 



, x r i d y\ 1 x r i d y l d y 3 , 

+ ^[z]-^ + ^6g id [z] — — + .. 
z = {x a ,v ai v a ,...},t,y; a = 1, 2, . . . , n, 



since it now moves along a geodesic on the spacetime 
metric generated by the rest of the n masses. There- 
fore, if L n +i is the (n + l)-body lagrangian less the 
M e (-1 + (1/2) (dy/dt) ), the deviation of the space- 
time metric from Minkowski 8g^ v can be read off the 
action of the test particle using the prescription: 



'l(7f)0 [t; •?] 



d 



dM, 



L n+ i[y = x] 



Sgoi [t, x 
Sgij[t, x 



d dL n+1 
dM e d(d y ydt) 

_d d 2 L n+1 

dM e d{dy*/dt)d(dyi/dt) 



y=y=---=M e =0 

[y = x] 



y=y=---=M e =0 



[y = x] 



y=y=---=M e =0 



We see that understanding and testing the dynamics 
- the equations that govern the time evolution of the 
{x a } — is intimately tied to understanding and testing 
the spacetime geometry g^ of the solar system. 

The outline of the paper is as follows. In section ILTl 
we set up a lagrangian description of the system of n 
compact astrophysical objects by approximating them 
as n point particles. Einstein's equations can then be 
solved pcrturbativcly as a Born series, whose graphical 
representation are the Feynman diagrams containing 
no graviton loops; the result of summing the diagrams 
yield the n-body effective action we seek. (Our descrip- 
tion will be brief because it will merely be an overview 
of the methods developed in [H[.) We then sketch the 
algorithm that could be used to generate the necessary 



gravitational waves produced by and interacting with the n 
compact objects. In this present paper, we shall focus only on 
the conservative part of their dynamics up to 2 PN. 
5 This observation can be found, for example, in Damour and 
Esposito-Farese 0], 



We work in cartesian coordinates and employ the ry^,, = 
diag[l, —1, . . . , —1] sign convention. The Einstein summation 
convention is adopted. Greek letters run from to d— 1 while 
English alphabets run from 1 to d — 1. 
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Feynman diagrams contributing to the effective action 
up to an arbitrary PN order. In section HTT1 we calcu- 
late the individual diagrams that occur at the Newto- 
nian thru 2 PN order and present the effective action 
of the n-body system up to certain integrals (f2"Tl l2"3l 
141) . As a by-product, we reproduce the known 2 PN 
2 body lagrangian. In the appendixes, we discuss the 
integrals encountered in the diagrams; the algorithm 
for generating the N > 2 graviton Feynman rules on 
a computer; and also display the Feynman diagrams 
occurring at the 3 PN order. 

II. THE n-BODY SYSTEM 

The assumption that we have a system of n com- 
pact objects, with their typical size r s much smaller 
than their typical separation distance r, i.e. r s <C r, 
suggests that the detailed structure of these objects 
ought not affect their gravitational dynamics, at least 
to leading order. These n objects could then be viewed 
as point particles. Because the most general action for 
a point particle must be some scalar functional of its 
d- velocity = dx£/ds a and geometric tensors 7 (and 
possibly the electromagnetic tensor F^, if large scale 
magnetic fields are present) integrated over the world 
line of the said particle; part of the action is already 
fixed to be of the form 



also ignore electromagnetic interactions.) Here, the 
{M a } lend themselves to a natural interpretation as 
the masses of the astrophysical objects and — M a J ds a 
describes a structure-less, mathematical point particle. 
At the 5 PN order and beyond, one would be com- 
pelled to include as many of the non-minimal terms 
as is required to maintain the consistency of the field 
theory up to a given level of accuracy. Physically, this 
means one has to begin accounting for the fact that, 
even if one neglects their rotation, astrophysical ob- 
jects are not really point particles and their individ- 
ual mass distributions and sub-structures do produce 
gravitational effects. To give the coefficients {c^} of 
these non-minimal terms physical meaning, one would 
have to compute (in-principle) measure-able quantities 
both in the actual physical setup and with the point 

particle terms in fl). The {<4?' > } are then fixed by re- 
quiring the results of the latter match the former: for 
instance, if we have multiple non-rotating black holes 
bound by their mutual gravity, then one could calcu- 
late the partial wave amplitudes of gravitational waves 
scattering off the Schwarzschild metric and match the 
point particle computation onto it by tuning the coef- 
ficients {M a ,c^} appropriately. 

Up to 2 PN, the gravitational dynamics of the n- 
body system, with x% denoting the /ith component of 
the coordinate vector of the ath point particle, is there- 
fore encoded in the action 5, where 



M a f ds a fl + c[ a) R^ a pW 

a=l J \ 

+ ci a) R llva0 R\ a p u v y a « 



+ 

CpR 



va.fi 



(i) 



where ds a is the infinitesimal proper time of the ath 
point particle and the ". . . " means one really has an 
infinite number of terms to consider, since the only 
constraints at this point are that each of them is a co- 
ordinate scalar and that none of them can be removed 
by a re-definition of either the metric g^ v or the photon 
field A^. 

However, as argued in [l3| . unless the n objects 
are very large or have very large dipole and higher 
mass moments, it is expected that the minimal terms 
{— M a J ds a } would suffice up to 4 PN order (see also 
§1.2 of Blanchet [|| for a discussion), and in what 
follows we will compute with them only. (We will 



S = Sqr + S pp 
Sgr 



Spp — 
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Moreover, we expect the metric of spacetime to de- 
part markedly from Minkowski only close to one of 
these n compact objects, where it is irrelevant for the 
problem at hand, and thus we can expand the metric 
about jy^: 8 



9nv Vfj>v ~l~ 



M. 



(d/2) 
pl 



7 The conventions for the Christoffel symbols T^^n, Riemann 



tensor R^ap = g^xR ^ a p: Ricci tensor and Ricci scalar 
1Z can be inferred from the formulae in appendix lAl 



The factor of M 



l-(d/2) 
pl 



ensures that the graviton kinetic term 



does not contain M p \. Also, for the rest of this paper, we will 
raise and lower indices with r)^. 
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The general relativistic effective lagrangian 
L e s[{x a ,Va,v a , . . .}} for n objects can now be 
computed via the prescription usually associated with 
perturbative quantum field theory, namely, as the sum 
of fully connected diagrams: 



exp 



i dt L e ff 



d-l ,. 

I Vh^ exp [iS + iSgf] 



exp 



S IT 



E /Fully connected\ 
\ diagrams / 



x ( d*h\p - -dfjh ) , 



h = /V\ 



(7) 



S G K[h 2 ] + S gi [h 2 } 



(9) 
(10) 



This choice corresponds to the linearized de Dondcr 
gauge j|"tW = d»h m - ±d a h^ = 0. 

The subscript "cl" (short for "classical") in {5} in- 
dicates the Feynman diagrams with graviton loops are 
excluded. As already remarked, evaluating these clas- 
sical Feynman graphs amounts to solving Einstein's 
equations for /i M „ via an iterative Born series expan- 
sion. 



A. Physical Scales In The n-Body Problem 



or alternatively, 



exp 



dt L 



eff 



: Af exp 
x exp 



iSi 



1 6 

i 8J^ 



,7=0 



-\fd d xj d d y jfDl . XT [x-y]jf 



(8) 



with D a p. XT [x — y] being the Feynman graviton 
Green's function and 5i[(l/z)(5/<5 J M "] indicates we are 
replacing every graviton field in ©, less the graviton 
kinetic term, with the corresponding functional deriva- 
tive with respect to J^ v with the same indices. In par- 
ticular, because the graviton field is symmetric in its 
indices, we have 



The expression in ([5]) is the functional integral ver- 
sion of the statement that, up to an (for current pur- 
poses) irrelevant factor AT, to compute the effective ac- 
tion, one needs to expand exp[iS'i] and, for each term 
in the series, consider all possible Wick contractions 
between the graviton fields. In the next subsection, we 
will use it as a guide to devise an algorithm for gener- 
ating the necessary Feynman diagrams at a given PN 
order. 

A gauge fixing term S g i (O has been added to make 
invertible the graviton kinetic term in the Einstcin- 
Hilbcrt action ([3]), whose explicit form then reads 



It is possible to begin computing the diagrams in © 
after only expanding @ in powers of graviton fields, 
performing a non-relativistic expansion afterwards and 
keeping the terms needed up to a given PN order. 
However, we will now show that it is more efficient 
if one also expands the action in terms of the num- 
ber of time derivatives and powers of velocities {v a } 
they contain, before any diagrams are drawn and cal- 
culated, as this will allow one to keep only the nec- 
essary terms in ([2]) such that every Feynman diagram 
generated from them scales exactly as v 2 ® , for a given 
PN order QeZ+. 

To this end, we note that, because we are assuming 
that the n objects are moving non-relativistically, with 
their typical speed v -C 1, we already know that the 
lowest order effective action must give us Newtonian 
gravity: 



S e S = 



(dt( £ \M a vl 
J ^ l<a<n " 



Rab 



l<a,b<r 



d-3 
ah 



with k being some (presently unimportant) dimcn- 
sionless number. This prompts us to associate with 
this lowest order action S c whenever this particular 
product of masses M, separation distances r and time 
occur in the action: 



dt Mv 2 



dt M 2 ~ d M 2 r 3 ' 



(11) 
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We may then obtain from pip 



S r ~ Mvr, 



r d—3 



(12) 



where the first relation holds because the only phys- 
ical length and time scales in the problem for a fixed 
coordinate frame are the typical separation distance 
r and orbital period r/v. Similarly, we relate all time 
and space derivatives and integrals to appropriate pow- 
ers of r and frequency v/r, 



dV 



(13) 
(14) 



where the 5[x° — x'°] relation will be needed shortly. 

Next, we observe that the real part of the Feynman 
graviton Green's function obtained from inverting ((9|) 
and (|10|) can be expressed as an infinite series in time 
derivatives: 9 



Re (0 | T { V [x° , x] Kp [x'° , £] } | 0) 

p °o F [ rf-3-2m | 

2-*> M 



m— 
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4 1+m 7rVrfl + m] 



d d 



\x - cc'\d-3-2m ^ dec' dx° 

oo 

E 



S[x° - x'°) 



iP^p ^ f d d -V e^-(g-g) 

/ d d 



Vdx'° dx° 



d-2 



to half its retarded plus half its advanced counterpart; 
see Poisson [ijj for a discussion. The introduction of 
an additional background field in [l3| takes into 
account the imaginary part of D^ v . a p[x — x'], which 
describes the dissipative part of the dynamics - the 
interaction of gravitational waves produced by and in- 
teracting with the n point masses. In this paper, we 
are focusing only on the conservative part of the dy- 



namics, and hence will ignore Im D a p\x — x'\. 

The zeroth order term in Re D^ v . a p[x—x'] ([T5j). with 
no time derivatives, can be obtained by inverting ()10D . 
i.e. the graviton kinetic term with only spatial deriva- 
tives. Diagrams involving the higher order terms in 
(fT5")) . with time derivatives, can be gotten by treating 
([5]), the graviton kinetic term with only time deriva- 
tives, as a perturbation. For instance, the first correc- 
tion to the Newtonian gravitational potential due to 
the finite speed of graviton propagation is proportional 
to / dtdt a dt b 6[t - t a ]S[t - t b }(d/dt a )(d/dt b )\x a [t a ] - 
Xb[tb]\ 5 ~ d , which could also be viewed as a contraction 
between two distinct world line operators of the form 

oo [x a ], from |g]), with one in- 
sertion of 

Keeping in mind that each diagram is built 
out of contracting graviton fields {h^ u [x\h a fj[x']) = 
(15) D^ u . a p[x — x'] from distinct terms in @, this im- 
plies every graviton field in (J^J) should be assigned a 
scale that is square root that of the lowest order non- 
relativistic Green's function containing no time deriva- 



tives. The \x- 



#13- 



„3-d 



dependence implies that 



spatial derivatives on h^ iv ought to scale with one less 
one power of r of the same. By recalling (|14[) , we then 
have 



d h^[x°,x] ~ r- d /V/ 2 



(16) 
(17) 
(18) 



where our notation alludes to the fact that the clas- 
sical Feynman graviton Green's function is the non- 
interacting-vacuum expectation value of the time or- 
dered product of two graviton fields. 

That there are only even number of time derivatives 
reflects the relationship that the real part of the Feyn- 
man Green's function for a massless graviton is equal 



Putting the scaling relations from (fTTj) , (jl2[) . (fT~3|) ■ 
(Til)) . (Til)]) , |[T7) and (JTHJ) into the action ©, wc then see 
that, upon expanding |(5J) in powers of graviton fields, 
velocities v a , and the number of time derivatives on 
hptv, each term in the action now scales homogeneously 
with S c and v: 



The momentum space representation in the following expres- 
sions — which is related to its position space counterpart via 
I IB31 - will be useful when contracting graviton vertices com- 
ing from the cubic and higher in terms in the Einstein- 
Hilbcrt action (J3J , because manipulation of spatial deriva- 
tives on hpiu become algebraic manipulation of momentum 
dot products in the numerator. 



J dt a O& [n, a, e s ]uf ~ si S v ^ + 2a+^ 

J d d x O v [m,tp] ~ 5 C 1_ f v 2m ~ 4+ ^ (19) 

where offl [n, a, es] denotes the world line term in 
dH) associated with the ath particle containing exactly 
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n graviton fields, less the v 1 " term from ^^v^v^, with 
a total of es spatial indices (for example, a term with 
hijtfv 3 hokV has = 3). Note that there is usually 
more than one term for a given es, so one has to sum 
over all possible terms. O v [m,ip] denotes the term in 
([3]) containing exactly m graviton fields (to > 2), with 
precisely ip time derivatives (ip = 0, 1, or 2). 

Given these results in (fH))) . and given rii v ) number 
of graviton vertices from ([3]), nr w ) number of world 
line operators from (j?]), and TV total number of gravi- 
ton fields (so that N/2 is really the number of Green's 
functions in the diagram) one can work out that every 
Feynman diagram in the theory arising from products 
of these operators must scale as 

so that such a diagram contributes to the {n^ w ) — 2 + 
As/2) PN effective action; and 7i(„) + — y = 1; 
as the non-relativistic expansion, for the conservative 
part of the dynamics, is a series of the schematic form 
S e g = Sq + S2V 2 + S^v 4 + . . . , where each term of the 
effective action has to contain the appropriate prod- 
ucts of masses, velocities and time integrals such that 
S n ~ S c , with the v n factored out. Here, As is a 
positive integer that is the result of summing powers 
of speeds coming from time derivatives, velocities con- 
tracted with graviton fields (such as hijV l v J ), num- 
ber of graviton kinetic terms with time derivatives 
© inserted, and the factors of arising from the 
term "q^v^v^ inside the proper time ds a - We have 
used, in deriving the exponent of v, the constraint 
that nr v ) + nr w ) — y = 1. The fact that no diagram 
can scale grea ter than the first power of <S C has been 
proven in [131 ]. Observe that, with only the minimal 
— ^2 a M a J ds a terms included, these scaling relations 
are independent of the number of spacetime dimen- 
sions. 



B. Algorithm 

We are now in a position to describe an algorithm 
that could, with an efficient implementation and suf- 
ficient computing resources, generate the necessary 
Feynman diagrams, for a given subset of world line 
terms in ([T]), up to an arbitrary order in the non- 
relativistic PN expansion. 

For a desired scaling (|20[) . corresponding to a spe- 
cific PN order, one can insert in $2$ explicit factors of 
iS c and v according to the results in (fl"9]l . so that one 
may employ Mathematica 10 [2(j to extract the relevant 



products of operators, i.e. pre-contraction, in the tay- 
lor series expansion of the exponential in the path inte- 
gral ([8]) using either its Series or Coefficient command. 
Observe that every fully connected diagram T>-p con- 
structed out of each term in the taylor series expansion 
of expfiSi] in (8J is the first non-trivial term of the se- 
ries expansion of cxp[2?p]- Hence the gravitational ef- 
fective action is the sum of all fully connected diagrams 
constructed from the terms in the series expansion of 
expfiSi]. Moreover, one does not need to include in 
the code the explicit form of the graviton fields, veloc- 
ities, etc., but it suffices to have placeholders, such as 
the ones used in (|19p . containing enough information 
to re-construct at the end the relevant types of gravi- 
ton fields considered (/loo, hoi or hij), which point mass 
the field(s) belongs to, factors of velocities, number of 
time derivatives in the m-graviton term(s), numerical 
constants from taylor expanding the exponential and 
the square root in the infinitesimal proper time, and 
so on. 

Next, the combinatorics of contraction can be han- 
dled by Mathematica, by assigning to every graviton 
field in a given product a distinct number, so that such 
a product corresponds to some list, say {l,2,...,s} 
for a product of s graviton fields. A permutation of 
{1, . . . , s} is equivalent to a Feynman diagram if and 
only if it leaves no numbers fixed and the resulting 
permutation operation 7? can be factored into products 
of disjoint 2-cycles, i.e. tt = (ai a2){a^ 04) . . . (a s _i a s ) 
with {ai, a2, . . • , a s _i, a s } being a re-arrangement 
of the original set {l,...,s}. For instance, the 
set {3,4,1,2} means one would have to contract 
graviton field "1" with graviton field "3"; graviton 
field "2" with graviton field "4" . (The requirement 
that each fully connected diagram scales as S^. en- 
sures there will be no quantum corrections.) What 
remains is removing those diagrams that arc not 
fully connected. One method of achieving this is 
to check if the permuted set {ai, a%, . . . , a s } can 
be factorized into two or more disjoint sets, where 
each of these individual sets contain only terms 
that are contracted amongst themselves. The term 
/ dt a J dt b (hoo[x a }h 00 [x b }) J dt c J dt e (h o[x c }h 00 [x e }}, 
for example, can be represented as {2,1,4,3} = 
{2, 1}{4,3}, i.e. factorize-able; whereas 

J dt a J dt b J dt c (h o[xa}h o[x b }}{ho [xa]h 00 [x c }) is 
equivalent to {3,4,1,2} and not factorize-able. Such 
a prescription can be implemented with a suitable 
adaptation of the command Permutations. 

Once the contractions are determined for a given 
product of terms from expanding cxpfiSi] , if the partic- 



We frame this discussion around Mathematica, but this al- 



gorithm can most likely be implemented with any software 
with similar symbolic differentiation and combinatorial capa- 
bilities. 
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ular diagram does not involve terms from the Einstein- 
Hilbcrt action, it can be computed automatically be- 
cause it would be built out of products of the gravi- 
ton Green's function. The scalar portion is \x a — 
xt\ 2s+3 ~ d , where s is the number of © inserted; while 
there will also be factors of velocities proportional to 

iP^v-aj3Va v a v b v b f rom P^l> - When insertions of ^ are 
present, one would have to take the appropriate time 
derivatives afterwards. (Some care needs to be exer- 
cised in keeping track of the time 5-functions when 
doing so - see ([2~4"]) for an example.) For diagrams 
with graviton vertices, although they may not be cal- 
culated automatically, the required Wick contractions 
and permutations of particle labels can be displayed so 
that the user does not have to figure out the combi- 
natorics manually, but rather focus only on the tensor 
contractions of the graviton Fcynman rules, manipula- 
tion of the momentum dot products and the ensuing 
Feynman integrals. Furthermore, some of the higher 
PN diagrams involving graviton vertices will be prod- 
ucts of lower PN graviton vertex diagrams with other 
expressions that can also be automatically calculated 
- such as factors of the graviton Green's function 
with or without insertions of © contracted into veloci- 
ties, (h^[xa]h a p[xb])v^v^v^v^ , etc. (See, for instance, 
Fig. ([5] a) and Fig.((5]c).) This implies the evaluation of 
such higher PN diagrams with vertices may most likely 
be automated if a repository of these lower PN gravi- 
ton vertex diagrams is kept. 

As an illustration of the utility of such an algorith- 
mic approach, we have generated in appendix [C] the 3 
PN Feynman diagrams for the minimal point particle 
action — Y] M a J ds a . We also maintain a web page 
at [23||, where the Mathematica code used in this paper 
can be found. 



C. n-body diagrams and superposition 

Now suppose one wants to calculate the lagrangian 
for n point particles up to the Qth post-Newtonian 
order. Then the exponent of v in ([20]) tells us that the 
maximum number of distinct particles that can appear 
in a given Feynman diagram is 

max[n( w )] = Q + 2 

and so at a fixed post-Newtonian order Q, obtaining 
the Fcynman diagrams for the (Q + 2)-body problem is 
sufficient for obtaining the lagrangian for the arbitrary 
n-body problem. In particular, at the 2 PN order, we 
see that the n-body problem is equivalent to the 4- 
body problem. For a general Q PN order, the diagrams 
for n > Q + 2 point particles can be obtained by sum- 
ming the diagrams for the n = Q + 2 case over all the 



particles in the system, since no additional distinct di- 
agrams are needed. For n < Q + 2 point particles, the 
relevant diagrams can be gotten from the n = Q + 2 di- 
agrams by setting the masses Mq +2 , Mq + \, . . . , M n+1 
to zero. Even with the non-minimal terms beyond the 
— ^ a M a J ds a included, it is apparent that superpo- 
sition will continue to hold at any given PN order once 
we have computed the effective lagrangian for a suf- 
ficient number of distinct point particles, since each 
Feynman diagram can only contain a finite number of 
world line operators. 



III. RESULTS 

We now present the diagram-by-diagram results for 
the computation of the effective action up to 2 PN. Be- 
cause the calculation is long and saturated with tech- 
nicalities, the reader only interested in the final results 
may simply refer to ([21]) for O[v ], {23]) for 0[v 2 ], and 
gl] S3 S3] [44]) for the 0[v 4 } effective lagrangians. 

All Feynman diagram integrals are evaluated within 
the framework of dimensional rcgularization, where the 
number of space-time dimensions, d = m — 2e, is some 
infinitesimal deviation from a positive integer: i.e., 
m € Z + and e = + . For some of the more difficult 
integrals encountered at the 2 PN level, we will restrict 
our interest to that of the physically relevant case when 
d = 4 — 2e. Within dimensional rcgularization, inte- 
grals such as / d d_1 p (p 2 )" CT , / d d_1 p p i (f)~ a and 
J d d ~ 1 p p l pi (p 1 )~ (J are set to zero. This can be justi- 
fied formally by setting to zero the appropriate a, p, 
or t exponent of ([B4]) . |B5]) . ([B6]) or (|B7]) . since T[z] 
diverges as z — ► 0. 

Because it is easier to manipulate momentum dot 
products than derivatives, both the Fcynman rules for 
the graviton vertices are derived and the tensor con- 
tractions of graviton vertices are performed in fourier 
space. (See appendix [A] for an algorithm that could 
generate the A-graviton Feynman rule for N > 2.) We 
will thus present the master integrals for each diagram 
first in momentum space. 

Notation A few words about the notation used: 
the time argument of the ath particle is t a , so that 
x a = a^a[^a]- However, if the spatial coordinate vec- 
tors {x a \a = l,2,...,n} and their time derivatives 
occur within a single time integral J dt, then it is 
implied that they all share the same time argument 
t. R a b = x a — xt and its Euclidean length is R a b = 
\x a - x b \ = - xl)(xi - x J b )) 1/2 . The partial 

derivative df = d/dx z a refers to the derivative with re- 
spect to the zth component of the spatial coordinate 
vector of the ath particle. The spatial velocity of the 
ath particle is v a = v a {t a ] = dx a /dt a = x a , and its 
acceleration is v a = v a [t a ] = d 2 x Q /d^ = x a . When- 
ever we compute in fourier space, the relevant sign and 
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FIG. 1: Newtonian gravity. 



7r conventions are encoded in the following definition: 
f[x] = (27r)~ d J d d p f[p] exp[ipo^°] cxp[— ip ■ x\, where 
/ is some arbitrary function, and x and p are its coor- 
dinate and momentum space arguments respectively. 

Feynman diagrams A blob with some letter "a" 
at its center represent a world line operator from (0J 
belonging to the ath particle, with the indices of its 
various graviton fields € {hoo, hoi, hij } indicated 
on the side, {a, b, c, e} are distinct labels. A line rep- 
resents the lowest order graviton Green's function with 
no time derivatives. The x on a line represent an in- 
sertion of (J5J). A black dot with k lines attached to 
it is the fc-graviton piece of (J3j) with zero time deriva- 
tives. The fc-graviton piece of ((3]) with 1 or 2 time 
derivatives will be indicated with a "1" or "2" respec- 
tively; see for example Fig.(j6]a) and Fig.([7|a). The 
v k appearing alongside the graviton indices of the ath 



world line operator indicates which power of u 2 from 
expanding — M a J dx°(l — v^ + ... ) 1 / 2 needs to be in- 
cluded. Every Feynman diagram displayed serves dual 
purposes: it represents the class of diagrams that can 
be obtained from it by permuting particle labels; but 
the result of the diagram shown in the body of the text 
is always for the specific choice of labels in the figure. 
(The exceptions are the diagrams where 2 3-graviton 
vertices are contracted: Fig.([S|c), Fig.flSjd), Fig. (fT2"] d) 
and Fig. (fT4] cV We will discuss the notations there.) 



(a) No permutations necessary. 

r\ 00 'irz\ 

(b) 2 permutations from swapping a <-> b. 
^00 00^ 

(c) No permutations necessary. 
^00 00^ 

0-. .0 



■00 

® 

(d) 2 permutations from swapping a <-> b. 

oo^oo oo ^ 

© © - © 

(e) 2 permutations from swapping o <-> b. 

A 00, v 2 00^ 

(f) 2 permutations from swapping a <-> b. 



FIG. 2: 1 PN two body diagrams 



present the results and not discuss any of the calcula- 
tion in detail. 



1. 2 body diagrams 
The 2 body diagrams are displayed in Fig.©. 



PN 



At the lowest order, we have Newtonian gravity com- 
ing from the single diagram in Fig.([T]) and the usual 
kinetic energy. 



Ka<n 



2 ^-8 r[ ^i] G ^ Ma M b 



- y 

a^b 



B. 1 PN 



R 



d-3 



(21) 



At 1 PN order, we have 2- and 3-body diagrams. 
Since the lagrangian at this order has been com- 
puted numerous times in the literature, we will merely 
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(a) No permutations necessary. 
^-sOO 00^00 00^ 

(b) 3 permutations: a, b, or c for the middle label. 
FIG. 3: 1 PN three body diagrams 

2. 3 body diagrams 
The 3 body diagrams are found in Fig. ([3]). 



Fig. (0a) = -i / dt 



(d_ 3 )2 r 2 [ ^3 ] MaMbMc 

128(d~2) 2 TT d - 1 M^ d ~ 2) 
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Fig. (0b) = i dt 
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R: 
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E 
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(23) 



Setting d = 4 recovers the known result in the lit- 
erature; for instance, equation (38c) of Damour and 
Schafer 0]. The d > 4, 2 body version of ^ has 
been computed by Cardoso et. al. (l5| . 



For later use, we note that the master integral for 
the 3-graviton diagram is 



2 PN 



l 000000 



[q,r, 



d-3 



M (d-2)/2 \d-2 
1V1 p\ 



n 



d^m\ vi+vi+Pi 

(2nY 



-Q ~9 ~9 

vivivi 



x exp [ipi ■ x q \t q ] + ip 2 ■ x r [t r ] + ip 3 ■ x s [t s ]} 
x(2TT) d - 1 S^[p 1 +p 2 +p 3 ] (22) 



At 2 PN, we have 2-, 3- and 4-body diagrams. We 
shall classify the diagrams according to whether they 
involve terms from the Einstein-Hilbert action, i.e. di- 
agrams with or without graviton vertices. Whenever 
there are time derivatives acting on <5-functions, for ex- 
ample (d/dt a )8[t a — i\, it is implied that integration by 
parts is to be carried out. To save space, we will not 
display the explicit result of differentiation. 



1. 2 body diagrams 



3. 0[v 2 ] Effective Lagrangian 



No graviton vertices The diagrams that do not 
involve graviton vertices are: 



Summing the first order relativistic correction to 
kinetic energy from the rj^v^v 1 ' in the infinitesimal 
proper time ds and the diagrams from Fig.@ and 
Fig-© hands us the 1 PN order, d > 4 dimensional 
ro-body effective lagrangian: 

r (l pn) 



Fig.|g|a) = i dt dt a / dt 



(d-3)r [^] 



5127r^-(d-2) 



M n M h 



d d 



M^ 2 R d ab 7 Vdiadtb 



5[t - t a ]S[t - t b ] 



Fig.flUb) =* / dt 



r [g] M a M b 

16tt 2 M pl R ab 
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s * 00 00^ 
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(a) No permutations necessary. 

a>' % 

(b) No permutations necessary. 

© ® 

(c) 2 permutations from swapping a <-» b. 

^0i,v 2 oi ^ 

(d) 2 permutations from swapping a <-> b. 

© ® 

(e) 2 permutations from swapping a <-> b. 

00, v 4 oo 
© ® 

(f) 2 permutations from swapping a <-> 6. 

© ® 

(g) No permutations necessary. 

© X ' (a) 

(h) 2 permutations from swapping a <-> i>. 

© x © 

(i) No permutations necessary. 
s-^00, v 2 w 00^ 

(j) 2 permutations from swapping a <-> 6. 

FIG. 4: 2 PN two body diagrams with no graviton vertices: 
1 of 2 
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(h) No permutations necessary. 

FIG. 5: 2 PN two body diagrams with no graviton vertices: 
2 of 2 
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Fig. (0c) 
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(a) 2 permutations from swapping a «-» 6. 
^-x 0/ 00^-x 



Fig-((5]f) contains a first order relativistic correction 
to the graviton Green's function. Integrating over the 
time 8- function with no time derivatives acting on it, 
before integrating by parts, the resulting integral be- 
comes 
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with a common time argument f for both x a and 
x b in the factor | . . . | 3 ~ d . 
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Graviton vertices The rest of the 2 body diagrams 
contain terms from the Einstein-Hilbert action. 

Note that the form of the fourier space master in- 
tegrals associated with each class of diagrams usually 
comes about after some manipulation of momentum 
dot products, application of the identity 2p a ■ p b = 
T(Pa T Pb) 2 ^Va^vi an d its analogs, and the use of 
momentum conservation X) r=1 p r = 0, i = 3 or 4. 

Fig.([6|a, b) The 2- and 3-body version of 
Fig.([6|a, b) requires the following master integral: 



-foioooo [q, r, s] 
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(d) 2 permutations from swapping a «-> b. 

FIG. 6: 2 PN 2 body diagrams with graviton vertices: 1 of 
4 
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PiPm 



X (2TT) d - 1 S d - 1 [p 1 +p 2 +p 3 ] 



(25) 



Notice that ap 3 s (with s = 1, 2 or 3) in the numerator 
can be obtained by differentiating the appropriate ex- 
ponential, i.e. p J s exp[ip s -x r ] = —i(d/dxl) exp[ip s - x r \. 
Our approach to the fourier integrals arising from this 
and the rest of the diagrams with graviton vertices is 
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(d) 2 permutations from swapping a «-> b. 
FIG. 7: 2 PN 2 body diagrams with graviton vertices: 2 of 
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(d) 2 permutations from swapping a <-> b. 
FIG. 8: 2 PN 2 body diagrams with graviton vertices: 3 of 



to first substitute the momentum (J-function(s) with its 
(their) integral representation(s), 



(27T) 



d-1 xd-1 



j d— 1 

d z exp 



-J :- 



(26) 



and next use (|B3[) to re-express the original mo- 
mentum integrals as position space ones, with the mo- 
mentum dot products in the numerator converted into 
derivatives on the resulting integrand. 

In this regard, the 2-distinct particles case usually 
requires more care than the 3- and 4-distinct particles 
cases. We shall illustrate this with Fig.([6|a), where 
x q = x r = x a and x s = Xb- The time derivatives 
occurring in Fig.dHja) are 
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= v{[t a 
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dxi [ta] 
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d^ 


= <[t b ] 


d 

dxi[t b ] 


d 

d^ 


= 4[t c ] 


d 

dx{[tX 



with each partial derivative acting only on the ap- 
propriate \x a — z\ or \xb — z\ with the same time argu- 
ment as the velocity vector contracted into it. 

When d = 4 — 2e, we therefore have 



Fig. (5| a) 
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After differentiation, these integrals can then be 
evaluated using (|B7|) . This leads us to 
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A similar analysis for Fig.([(5]b), making use of (|B4 
and (|B6[) . gives 
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Before proceeding further, it is useful to introduce 
the following master integral that would occur in sev- 
eral 2-, 3- and 4-body Feynman integrals: 
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(d) 6 permutations. See text for discussion. 

FIG. 9: 2 PN 2 body diagrams with graviton vertices: 4 of 
4 
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(27) 



where we have provided both its fourier and posi- 
tion space representations. In appendix [Bl we obtain 
I 3 [a, b, c] in closed form when d = 4 — 2e, up to O[e ] 

(EH. 

Fig.([6|c,d) We now turn to Fig.([6|c) and 
Fig.([6|d). Its associated master integral is 
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d-3 (pV # g ) 2 + (p 2 ■ t7 g ) 2 4- (p 3 ■ 
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2(d-2) 

x exp [i (pi ■ Xq [t q ] + p 2 ■ x r [t r ] + pz ■ x s [t s ] )] 

x (2TT) d - 1 5 d - 1 [p 1 +p 2 +p 3 } (28) 



The term proportional to may be integrated 
in arbitrary d-dimensions by integrating over the 
momentum that is absent in the denominator (af- 
ter cancelation), followed by an application of 
()B3p . because it reduces to a product of the form 

Usi 2 *) 1 "*! d^psC^iPs ■ Rsu]/f s . This type of 
fourier integral will occur frequently. 

The second term containing momenta dotted into 
velocities has the position space representation 
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terms, one with a R l az R J az in the integrand and the 
other with S tJ . If we simply set D = d and employ 
the formulae (|B4j) and (|B6j) . each of the two terms will 
be ill-defined. However, displacing D = d + n in (|B4|) 
and (|B6[) . and performing a laurent expansion in k 
afterwards would yield a finite result that is not trace- 
less and is furthermore consistent with first doing the 
scalar integral (|B4j) with a\ = p and a 2 — d — 3 and 
then carrying out the double derivatives at the end. 

If one is interested only in the higher PN 2-body 
problem, the existence of such subtleties may be reason 
to stay within fourier space as far as possible for the 
evaluation of Fcynman integrals. For the 2 body case 
of the integrals (|29|) and analogous ones below, to avoid 
this subtlety for the terms where there are 2 factors of 
Raz d ( or 2 factors of R\ z d ) and the derivatives are 
acting on the R b i z d (or R^~ d ), we shall first do the 
integral using I3, before differentiation. 

Along this line, we further remark that 



d d - 1 zRH^d b J Rl: d 



can be integrated and then differentiated, whereas 



There is a subtlety when taking double spatial 
derivatives on a single factor of the Euclidean distance 
raised to the 3 — d power occurring within the Feyn- 
man integrals, such as dfdjR^~ d . Strictly speaking, 
because R^ LZ ~ d is the Green's function of the spatial 
laplacian operator S^d^d®, one needs to employ the 
formula 



d d 



dx i dxJ 47r( d - 1 )/ 2 |x| d - 3 



27r (d-l)/2 



(d-1) 



XX 



•|d+l 



where there is a <5-function term in addition to those 
following from straightforward differentiation so that 
one would obtain the correct result upon taking the 
trace of both sides. Insofar as the Fcynman integrals 
are concerned, however, it appears the J-function term 
may be dropped as long as proper regularization is 
used. For instance, if we try to compute the integral 



by first carrying out the differentiation (without in- 
cluding the S- function term), we would obtain two 



„ r)3 — d p3 — d r\b ob r>3 — d 
zH az Rbz d i d j H bz 



has to be differentiated first. One cannot begin with 
3 distinct coordinate vectors {x a , Xb, x c } in this inte- 
gral, engage ^3, and then set x c = Xb'- the last step in- 
volves terms like \im c ^b Rl c / Rbc and lim c ^&i?^ c 2 and 
hence is ill defined. 

We now employ ()B3|) . derivatives on I3, (|B4[) and 
(lB6l) to deduce 



Fig. (0c 



8M 



Ma Mb T , 

(3/2)(d-2) lijoooo[a,a,b\ 



pi 



d=4 



at 



M?M h 



1024M p > 2 ^ 6 



x (-3(iL • v a f - i(R ba ■ v a ) 2 + v 2 a i? 2 fc ) 



and 



Fig. (Eld) 



M n Ml 



(3/2) (d-2) AiOOOO 
Pi 

j. M a Ml 
2048M p V 2 < b 



16M 



d=4 . 

= I 



[a, b, b] 
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x {-4(R ab ■ v a f - 3(R ba ■ v a ) 2 + 5v 2 a R 2 ab } 



Fig.([Zla) It turns out Fig.(j7|a) and its 3-body 
counterpart are zero in d = 4 spacctimc dimensions. 
Because it is not apparent, we display the results here 
for the 3 body case. In terms of the master integral in 
(|2"T)l . we have 



Fig.ccgc) 

= -i J At J dt a J dt b J dt, 
(d-4){d-l)M a M b M c 



16(d - 2) 2 Mp! (d ~ 2) 



I 3 [a,b,c] 



d\ 2 /d\ 2 /d^ 2N 



,dt a J \dt b J \dt r 
x S[t a - t]5[t b - t]6[t c - t] 



(30) 



To be sure, when the number of distinct particles 
changes from 3 to 2, the integrals occurring in Fig. (J7|a) 
would be different from that in (|3"0"|) - which really is 
the result for Fig. (|ll| c) - but because it remains finite, 
Fig- (0a) s till vanishes due to the coefficient (d — 4). 

Fig.([7|b) Fig. (|7jb) corresponds to the first rcl- 
ativistic correction to each of the lowest order gravi- 
ton Green's functions in the 1 PN 3-graviton diagram 
Fig. (J2]d) . Its associated master integral is 



^oooooo x [q, r, s] 
I dt / dt r . 



dt r / dt 



(d-3) 2 



(d - 2) 2 M^~ 2)/2 



n 



d d -y 

d-1 



(27T) 



pj+pj+pj 

-A -O S) 

vivivi 



S[t - t r ]S[t - t s ] 



; S[t - t q 



dtp- 

j2 



^+^±^S[t~t q ]5[t-t s ] (^6[t-t r ; 
P1P2P3 \ dt i 

pj + pj+pj 

-9-9-4 

pivm 



8[t - t q ]5[t - t r ] 



M S[t-t s 



X exp [i (pi ■ Xq [t q ] +p 2 ■ x r [t r ] + p 3 ■ x s [t s 
x (27r) d - 1 5 d - 1 [p 1 +p 2 +p 3 ] 



(31) 



For the case of 2 distinct particles, ioooooox [a, a, b] 
becomes 



-foooooox [a,a,b] 



d=i 



dt I d 3 - 2 ^- 1+2£ i? b - 1+2£ 



128tt 3 M p1 

x (vid? + vi4d?d?) iC +2e 



dt / dt n / dt h / dt 



2M pl 

x 5[t - t a ]5[t - t b ]S[t - t c ] 

d 3 ~ 2e pi d 3 ~ 2e p 2 ( e^-^^-x^t.X) 

( 27r )2(3-2e) y 
iVg [tq] -Pl-Vg [tq] ■ Pi Vg [tq] ■ Pi ift .(S a [t a ]-x a [t b ]) 



I 



—1 
Pi 



647T 2 M p ii? ah 
2 



{v\\d? + v{vid? d!)R ab 



4M 



pi 



(vtd> + vivid»d»)l 3 [a >a ,b}. 



The first position space integral can be evaluated 
using (|B"4")) , (|B5)) and (|B6|) . The fourier integral after 
it vanishes upon integrating over the time (5-functions 
because the exponential becomes unity. We thus have 



Fig. (0b) 

M 2 M b 



16M 



\(d-2) 



^000000 x [a, a, b] 



pi 



2048M p V 2 i?^ 



x \ h{R ab ■ v a y + A{R ba ■ v b y 

- (4i?ab • 'v a + 2R ba ■ ij b + 3v 2 a + 2v 2 )R 2 ab 



Fig.([7|c,d) The master integral associated with 
Fig.(Hc) and Fig.([7|d) is 



loiojoo [q,r, s] 



I dt I dtq I dt r 



dt. 



5[t - tg}S[t - t r ]8[t - t s ] 
4(d-2)M^ 2)/2 



n 



A d-1 

d pi 



(2tt 



\d-l 



X < V ■ V r 



(d-4M+0l) + (d-2)0l 



-9-9-9 

P1P2P3 



•2(d-4)i 



-9-9-9 

P1P2P3 

^ P3 ■ Vq Pi ■ Vr + P2 ■ Vq P3 [ Vr 
T z l« -Jj -12-11-11 
P\P\P 3 

X CXp [i (pi ■ Xq[t q ] +p 2 ■ X r [t r ] + p 3 ■ X s [t s ])] 
X (2TT) d - 1 S d - 1 [pi+p 2 +p 3 ] 



(32) 
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A similar approach to the one taken for (|25| and 
([28]) . together with the integrals JB6| and JB7}, then 
yields 



Fig.(Hc) 

M 2 a M b 



x pl 



d=4 



df 



M 2 a M b 



512M p V 2 < b 



x ( v a ■ v b R 2 ab + R ab ■ v a R ab ■ v b - 4i? ha • v a R ba ■ v b 



dt 



MlM b 



512(d- 2)3My- 2) 7 rt( d - 1 ) J R^- 3) 



Fig. (Eld) 
df- 



r [^i] Ml Ml 



64(d - 2)M 



pi 



'000000 



[a, b, b] 



= -i / dt 



r[i 



Ml Ml 



256(d - 2fM^ d - 2 \W-^)R^-^ 



Fig.([9|a,b) The associated master integral for 
Fig. (0a) and Fig. (0b) is 



as well as 



d=4 



M?M b 



128M p Vi^ 



Fig.(g|a,b) Fig.|(l]a) and Fig.jS]b) are, up to nu- 
merical constants, products of v 2 with the 1 PN 3- 
graviton diagram Fig.((2]d), namely, 



Fig. ©a) 



16M 



[(d-2) 



Pi 



(If 



dt i%Ioooooo[a,b,a] 



(d - 3) 2 r [^-] MlM b ^ 
256(d - 2)2M p 2 f- V-ii^ d - 3) 



and 



Fig.©b) 



32M, 



[(d-2) 



pi 



df « 2 /oooooo[a, a, 6] 



(d-3) 2 r [^] 2 M 2 M b 



512(d- 2fM 2 J [ d - 2 \ d - 1 R 2 } b d ~ Z) " 



Fig.(g|c,d) Fig.(H]c) and Fig.flS]d) are, up to con- 
stant factors, the product of the lowest order graviton 
Green's function (/loo^oo) with the 1 PN 3-graviton 
diagram Fig.([2]d), namely, 



Fig.JSJc) 
= fdt- 



r [^i] MlM b 



128(d-2)M p f d - 2 V^i? d - 3 



'000000 



a, a, 61 



ah 



'00000000 



[q, r,s,u] 



I dt I dt q I dtr I dt„ 



dU 



x S[t - t q }5[t - t r }S[t - t s ]S[t - t u ] 



n 



d P£ 

d-1 



(27T) 



(d-3)(d(7d-51) + 86) 
X 4(d - 2) 3 M^f 2 

P1P2P3P4 

x cxp [i (pi ■ Z + p2 ■ x r + p s ■ x s + Pi ■ Xi 
x {2n) d - 1 5 d - 1 [p 1 +p 2 +p^+p i } 



(33) 



Employing (|B3|) hands us 



Fig. (0 a) 

MlM h 



96M 
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df 



(d - 3)(d(7d - 51) + 86)r 3 [**=2] M Q 3 M 6 
24576(d-2)3M p 3 1 (d - 2) 7rK £i - 1 )i?^ 3) 



and 



Fig.(0b) 



M 2 M fc 2 

a ^Z2T ^00000000 [«! a , 6] 



64M 



2(d-2) - 
Pi 



Fig.([9|c,d) The evaluation of Fig.(0c) and 
Fig.(0d), involving the contraction of 2 3-graviton ver- 
tices, requires the most effort. Its associated master in- 
tegrals in momentum space, after substantial algebraic 
manipulation reads: 
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/oooo-oooo [q, r,s,u} = /qooo-oooo [?> r > s > u ] + J oooo-oooo [l> r 7 s,u}+ /qooo-oooo [?> r > s > M ] 



^0000-0000 r : s i u ] 



= i dt dt g dt r / d£ s / dt u 



8(d-2) 3 M^" 2 



n 

^=1 



A d-l 
d-l 



(27T) 



-*2 ^9 ^9 -9 



(d-3) 2 (5d- 18) (pf +f 2 +f 3 +f 4 ) + (3d 3 - 40d 2 + 151d - 174)(pi+p 2 ) s 



(d-3) 2 



(4(d-2)(p>1 +pfpt) - (7d- 22)(pf +p1)(p1 +p\)) 



(pi + P2) 2 

x exp [ipi ■ x q [tq] + ip 2 ■ x r [t r ] + ip 3 • x s [t s ] + ipi ■ x u [t u ]] (2ir) d ~ 1 8 d ~ 1 
5 other permutations of {x qi x r , x s , x u } 



4 



'0000-0000 



[q, r, . 



dt / di, / dt r / dt s / dt u 



(d - 3) 2 6[t - - t r ]6[t - t s ]<5[£ - i„] 



4(d - 2) 2 M^f 2 



n 

^=1 



(2tt) 



d-l 



pi +p! +p1 +p! 



d-l rd-1 



P1P2P3PKP1 + P2) 2 

+ 5 other permutations of {x q , x r , x s , £«} 
-^0000-0000 [li r i s ; u ] 



exp [ipi ■ x q [t q ] + ip 2 ■ x r [t r ] + ip 3 ■ x s [t s ] + ip4 ■ x u [t u }} (2w) S 



' 4 

W—l 



dt I dt q I dt r I di„ I dt, 



(d - 3) 2 S[t - t q ]6[t - t r ]6[t - t s ]6[t - t u ] 
' 2(d - 2) 2 il^f 2 



n 

Ki=l 



(27T) 



d-l 



Pi -P3 P2- P4+P1- Pi P*2 ■ P3 



d-l rd-1 



pIpIpIpKpi +P2) 2 

5 other permutations of x r , x s , f u } 



exp [ipi • x q [tq] + ip 2 ■ x r [t r ] + ip 3 ■ x s [t s ] + ipi ■ x u [tu]] (2tt) 5' 



(34) 
(35) 



(36) 



(37) 



When contracting the 2 3-graviton Fcynman rules in 
momentum space leading to (|34I35I36I37|) . one makes 
a particular choice for the labels on the point particles' 
coordinate vectors. If the two particle labels on the 3- 
graviton vertex on the left hand side are {a, b} and the 
two on the right hand side are {c, e}, so that such a 
choice can be denoted as either (ab\ce), (ba\ce), (ab\ec) 
or (ba\ec) - there is a symmetry obeyed by the labels 
on either side - then the 6 permutations to be summed 
over in the definitions of /oooo-oooo. ^oooo-oooo and 
/qooo-oooo arc: 2 x (ab\ce), 2 x (ac\be) and 2 x (ae\bc). 
The particular permutation displayed in Fig.([9]c) is 
(aa\ba). It also represents the sum of the 6 permuta- 
tions, with each of the 6 terms giving the same result. 
In our notation, the sum is 6 x (aa\ab). As for the two 
diagrams in Fig.([9]d), they represent the sum of the 6 
permutations: 4 x (ab\ab) and 2 x (aa\bb). 

Turning to the integrals themselves, 
^0000-0000 [<7j r i s i u ] can t> e calculated in d-dimensions 



with (|B3[) . For the term containing (pi + P2) 2 
in the numerator, one would first sum over the 6 
permutations of the coordinate vectors {x q }; this is 
equivalent to holding fixed the {x q , x r , x s , x u } and 
permuting the names of the integration variables 
{pi, . . . ,Pi}. Upon doing so, one would find the rele- 
vant part of the integrand can be replaced as follows: 

(pi+P2) 2 /(p1»1) - 2(p1+p1+p1+p1)/(p!p1p1p1). 

After suitable re-definitions of the momen- 
tum variables of the form p^ = p a + Pb, the 
(P1P2 +Ptpt)/(fi +P2) 2 portion of /0000-0000 fe> r > s > u ] 
is, after integration, proportional to 5 d ~ 1 [x q — x r ] 
or 5 d_1 [0], i.e. zero within dimension regularization; 
while the (pf +pl)(p| +p1)/(pi +P2) 2 integrand again 
takes the (|B3|-form cxp[iJ2lp r ■ Rrul/iYllp 2 }- 

/qqq _qqqq[(7, r, s, u] can be dealt with by first inte- 
grating, for each of its four terms, over the momentum 
appearing in the numerator. We demonstrate this with 
the first term, dropping all numerical constants and 
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suppressing the time arguments. 



-4 

pi 



x exp • x q +p 2 ■ x r + p 3 ■ x s +P4 ■ Xu) 
x {2n) d ~ 1 5 d - 1 [p 1 +p 2 +p 3 +p A ] 




-0,-9-9 

.vinvi 



ip2 ■ (x r -S q )+W3 ■ {S s S q ) +ipi ■ (z„ — x q ) 



9A p 2 r iP2-(£ r -S a ) (P3+P4) 

u y -2 e — o — o / — * 1 — * \ 

p 2 Piplm+Pi) 



X e 



ip3-(Xs—£q)+ip4-(Xu.—£q) 



(38) 



have 2 distinct coordinate vectors occurring in each 
of them, i.e. (ab\ab). For the former, one may sim- 
ply use (jB4}, (Ell, (|B7]l . together with the cosine rule 

i?a,-iL =-(l/2)i^ b + (l/2)i?L + (l/2)i?L- For the 
latter (a6|a6) case when d = 4 — 2e, one can first inte- 
grate over y using ^3 , carry out the differentiation, and 
by applying the cosine rule the (R az + Rbz + Rab) ap- 
pearing in the denominator right after differentiation 
of J3 will be removed. The integrals that remain are 
tractable via (|B4|) . 

Summing up the contributions from the 6 permuta- 
tions for each of the Fig.([9]c) and Fig.([9]d) now gives 



192M 



2(d-2) 
Pi 



0000-0000 [a, a, a, 



d = 4 -i I dt 



MlM b 



16384M*tH»J&' 



where we have dropped the term proportional to 
/ & d l p2 exp[ipV (x r [t r ] — which is zero even if 

x r — x q , within dimensional rcgularization. The first 
piece in the second equality containing only squares of 
momenta in the denominator can be done with (|B3[) . 
For the second piece, the P2 integral is a {\/i)d\ on 
(IB3I) ; it is zero if x r = x q . The P3 and P4 integrals 
translate to an appropriate derivative on I 3 [a, a, b] or 
Iz\a,b,b] if x s = x u and x s 7^ x q . 11 If x s = x q and 
%u 7^ x*q ( or vice versa), the P3 (p^) integrals take the 
form of (|B5|) . and the remaining p 4 (p 3 ) integral is then 
a derivative on (|B3|) . The p 3 and p*4 integrals return 
zero if x 8 = x u = x q . 

What remains is ^0000— 0000 [*?' r ' s ' First replace 
the {pi + P2) 2 in the denominator with q and intro- 
duce a (27r) 1 - d J d d " 1 g^ d - 1 )[g -pi -p 2 ]. Applying 
(|2"rj)) on both the (5-functions then tells us that, in its 
position space representation - again ignoring the con- 
stant factors - -foooo-oooo [<!■> T -> s j u ] becomes 



/ d d_1 z 6 ij 6 mn 

- rf )(^^- d )(^- d )(5|ii s 3 - rf )(^^- rf ) 
{diR q y d ) {d r m Rl~ d ) (R^~ ) (d^R 3 s ~ ) {d™R 3 u ~ d )) 



x {(d!R qy 



For the 2 body problem, we would cither have one of 
the y or z integration involve only one of the coordinate 
vectors x a or x*b, i.e. (aa\ab) and (aa\bb) = (bb\aa), or 



11 The p3 and p4 integrals with the numerator p| + p\ removed 
are I^[q,s,u]s because one can recover the form H27H if one 
introduces an additional variable q = p% + p 4 and a corre- 
sponding integral and momentum conserving <5-function. 



and 



Fig.pd) 
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2(d-2) 



^0000-0000 a, b, b] 



d=4 



MlMl 



16384M p V 3 i?^ b 



2. 3 body diagrams 

No graviton vertices The 3 body diagrams that do 
not have graviton vertices are: 



Fig. (EH a) 



(d- 3) r[^r 

256(d-2) 2 7r d - 1 

MqMbMc 

M 2 /- 2) Rt 3 Rt 3 ' 



'be 



Fig.(dD|b) 



df 



(d-3)r [^] 2 

128(d-2) 7 r d - 1 
M a M b M c 



Fig.(JTU|c) = i dt 



(d-3)r 



'd-31 



128(d-2) 2 7r d - 1 

M a MbM c _ 2 

M3 d - 2) Ric 3 Rt 3Vc 



V 



Fig.dinid) 



df 



3r[¥] ; 



128 (d- 2) 2 tt 



2_d-l 
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© 



00 



00^.00 
-(c)--- 
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(a) 6 permutations: a, b, or c in the middle. 
Two choices for remaining two labels. 

^00 00^0/ Qj A 

® ' © • Cs) 

(b) 6 permutations: a, &, or c in the middle. 
Two choices for remaining two labels. 

00^ ij 00^, 

(c) 3 permutations: a, b, or c in the middle. 
00 00^-. 00, i/ 2 00^-. 

(d) 3 permutations: a, &, or c in the middle. 

A 00, v 2 00 00 00^ 

® " © • © 

(e) 6 permutations: a, 6, or c in the middle. 
Two choices for remaining two labels. 

^00 00^00 00^ 

(f) 3 permutations: a, b, or c in the middle. 

^00 00^ 

v00 
'000 
00; 

:00 

(fi) 

(g) 6 permutations: a, 6, or c in the center. 
Two choices for remaining two labels. 



00, 



00 



r~\00 



(h) 3 permutations: a, b, or c at the two ends. 



00 



r~\00 



r~\00 



® 



' 00^ 00^ 00 

(i) 6 permutations: a, b, or c for the repeated label. 
Two choices for remaining two labels. 

FIG. 10: 2 PN three body diagrams with no graviton ver- 
tices. 
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Fig. OH e) 
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128(d-2) 2 7r d - 1 
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As with its 2 body counterpart, Fig.fTUjf) requires 
some caution when taking the time derivatives. As an 
example, the integral in the a — x — c — b diagram, 
without the constant factors, is 



/d t „/da[« -«j 

'd\x c [t c ] - x b [t c }\ 3 ~ d d\x a [t a ] - f c [t c ]| 5 " d 



+ \x c [t c ] - %b[tc]\ 



dt c dt a 

3— d^ \ x a[ta] ^c[^c]|^ 



dt a dt c 



Fig.|T0]f) is the sum of a — x — c — b and b — x — 
c — a, but since differentiation and piecing together 
the relevant constants are straightforward, we will not 
display the result. 



Fig. p] g) = i dt 



1024(rf-2) 3 7rt(' i - 1 ) 
M 2 a M b M c 

M 3( rf -2 )i? 2( d -3)^-3 



Fig.fl^/d ^ Jj^r 



512(d-2)3 7 r3(' i - 1 ) 
M 2 a M b M c 

M3(d— 2) pd— 3 Tjd— 3 rid— 3 
pl ^ah H ™ K bc 



Fig.(TJJ|i) = « / dt 



512(d-2) 3 7r*( d - 1 ) 
M a 2 M fc M c 

j f3(d-2) R d-3 R 2(d-3) 
M pl K ab 



Graviton vertices The rest of the 3 body diagrams 
contain graviton vertices. 

Fig.(HI]a) Fig.frpa) requires ioioooo from (J23J). 
For the 3 body case, one sees that /choooo can be ex- 
pressed in terms of time- and space-derivatives on I3. 
Specifically, 



^oioooo [a, b, c] 



2(d - 2)M^f" 2)/2 



di / df Q / dt 6 / dt c 



x [* — t a ]*[t — * 6 ]*[* — * c ] 

x ({(2(d - 4)vid\ + (d- 5)v* a d?) h[a, b, c]} 



x — S[t b - t]S[t c - t] 
dt b 

+ { (2(d - 4)<Sf + (d - 5)<^) J 3 [a, 6, c]} 



x —6[tb-t]6[t c -t] ), 



which means 
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(a) 3 permutations: a, fo, or c carry the Oi indices. 
/~ \ 00 00/-^ 

©.. .© 



00 



© 



(b) No permutations necessary. 

/ — N 00 oo^ 

©-. -© 



©' 



(c) 3 permutations: a, b, or c carry the ij indices. 
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(a) No permutations necessary. 



© 



00 



00 



© 



+ (b—c) 



00 



G) 00 



© 



^ 00 v 

(b) 9 permutations. See text for discussion. 

®oo 00® 



Oi 



© 



(d) 3 permutations: a, b, or c carry the 00 indices, 
^-x 00 00 /~N 



(c) 3 permutations: a, b, or c for the repeated label. 

©00 oo© 



© 



00, © 



(e) 3 permutations: a, b, or c carry the v . 

FIG. 11: 2 PN three body diagrams with graviton vertices: 
1 of 2 
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Fig.([TT|b) As discussed for Fig.([7]a), we have a 
vanishing diagram when c? = 4 for 
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(d) 6 permutations. See text for discussion. 

FIG. 12: 2 PN three body diagrams with graviton vertices: 
2 of 2 



Fig. CU b) d = 4 
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Fig.([TT]c) The master integral for Fig. pTj c') is 
hjooao in (|28)) , whose 3 body expression in terms of I3 
is 



lijoooo [a, b, c] 



at 



16n d - 1 (d-2)ni^ 2)/2 



-J2 I r>3-d r/3-d 
V a{ H ab H ac 



(d-3)(Rl- d Rl d + R 3 c - d R^ d )) 
d-3 



2(d-2)M^ 2)/2 a a 



x (d?dj + d\d) + dldf) I 3 [a, b,c]\, 



so that 



__ . M a M b M c . , . 

Fig-dHc) = a^y%ooooKo,cJ 



8M 



pi 



f r>3 — d c>3— d 1 r>3— d t?3— tZ , r?3— d r>3— d\ 



Fig. ([12] a) The 3 body master integral for 
Fig. (HI] a) is /oooooox from (J3TJ) : 



-^oooooox [a, b, c] = I x [a, b, c] + I x [b, a, c] + I x [c, b, a] 

Ix [a, b, c] 

i(d-3) 2 



(d - 2) 2 M^~ 2)/2 
x [/ 3 [o,6,c] 



dt / dt„ / dt h / dt, 



r r^l r [^1 

647r rf-l { K ba R bc + K ca Rcb ) 



x ^6[t-t a ]6[t-t b ]6[t-t c 



In terms of /oooooox [a, 6, c], 



Fig. tlllj d^ The 3 body master integral for 
Fig. pTj d) is Ioiojoo from (p{2")) , which reads 



JoiOjoo [a, &j c] 

=*/<*( rm2 "",? 2)/ . 

x((d-4) (i?^- d + i?LX 3 c^) 

+ (d-2)i? c 3 - d i? c 3 ro 
1 < j 

2(d-2)M(^ 2)/2^;a?;^, 

x ((d-4)a?a; 



(d-3)(9^ + ^))/ 3 [a,6 jC ] 



leading us to 



MaMbMc T r , , 
Fig-CUd) = 2^y/oi0j0oKo,q 



2M 



pi 



Fig. (|ll| e) Fig. pTj e) involves products of z/ 2 with 
the 1 PN 3-graviton master integral /qooooo- 



Fig. (HH a) 
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Fig. (|12| b) For the class of diagrams represented 
by Fig.([T2|b). the 9 permutations include both types 
of diagrams where the 2 world line sources attached to 
the 3 graviton vertex not contracted with an additional 
world line can cither belong to the same particle (6 of 
them) or 2 distinct particles (3 of them). (Tig. ([12] b) 
itself belongs to the latter.) We observe, as we did 
in the 2 body case, that these diagrams are products 
of the lowest order (/loo^oo) with the 1 PN 3-graviton 
master integral /oooooo m (|22"|) . 



Fig.CTJb) 



r [jg] M 2 a M h M c 

Ud-Z) 



dt 



64(d- 2)tt— M p 3 j 
x (^ac ^000000 [a, b, a] + R^ b d ioooooo[a, c, a]) 
t (d -3) 2 r [i^] 2 T[^]MlM b M c 
1024A4 (d ~ 2) 7ri( d - 1 )(d-2) 3 



= -i dt 



( P 3 - d R 2 
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2(3-d)^3- 



(39) 



Fig. die) = 



MaMbMc-a T , , , 
dt ,,, „, «„/oooooo [a, 0, q 



16M 



(d-2) 



Pi 

, (d-3) 2 r[^] 2 M Q M b M c ^ 2 
-« / dt ±-±-i _ v 

256^- 1 (d-2) 2 M 2 1 (d_2) 



For the other 6 permutations where the 2 world line 
sources attached to the 3 graviton vertex not con- 
tracted with an additional world line belong to the 
same particle, the analog to the above 2 terms in (|39|) 
have different numerical factors in front of them. 
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Fig.([T2]c) Fig. p2| c) is a straightforward appli- 
cation of /qooooooo in (|33j) and (|B3|) . 
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^bc ^ab 
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Fig.([T2]d) Referring to the discussion under the 
2 body counterpart of Fig. (fT2] d), if a is the repeated 
label for a given 3 body diagram, then the 6 permu- 
tations of particle labels are: 2 x (aa\bc), 2 x (ab\ac), 
and 2 x (ac\ab). 

The master integrals for Fig. (fT2"] d') can be found 
in (|34I35I36I37[) . ^oooo-oooo ^ s the linear combination 
of products of (|B3|) . As we did in the 2 body case, 
Iqooo-oooo can be done using (IB3|) . its derivatives, and 
derivatives on I3. When a is the repeated particle la- 
bel, ^0000-0000 i s 



^0000-0000 [°i a i c ] 



d=4 



256M2j7t 3 



r2 _ e>2 p2 , n>2 _ p2 

^ I ^ab n ac ^ab ">" n ac ^bc 



-^ab + R lc - R lc , R ab ~ Rbc 



2R*m 



b L ac 



RficRbc 



Rac - Rbc Rlc ~ R\i 



RlbRbc 2Rl b Rl 



(40) 



We will leave -foooo-oooo[<7> <li T -> s ] f° r possible future 
work. 

In terms of the /0000-0000 m (|34I35I36I37|) . we have 
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64M, 
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3. 4 body diagrams 
No graviton vertices The vertex-less diagrams are: 
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(a) 4 permutations: a, 6, c, or e for the center label. 



(b) 6 permutations. See text for discussion. 
FIG. 13: 2 PN four body diagrams with no graviton vertices 
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™v-/ „o(e}-" -~ ™OJ + (a**b) 
00^^ 00 v ~ y 00 v - y 



For the class of diagrams in Fig. (fT3] b). there are 12 
ways to choose 2 out of the {a, 6, c, e} for the mid- 
dle two labels, but there is a reflection symmetry, and 
hence there are 6 distinct permutations of the particle 
labels. 



Fig.dElb) = i / dt 



r[*?] 



M a M b M c M e 



512(d- 2) 3 tt 
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3(d-2) 
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X i? 



3— d z?3—d Tj'3—d 

ae 



1 no- d po 

-"be ^ 



+ ^ac ^be ^ 



Graviton vertices The rest of the 4 body diagrams 
have graviton vertices. 

Fig.(fT4|a) The class of diagrams in Fig. ([T4"| a). 
like those of its 2- and 3-body counterparts, involve 
products of the lowest order (/loo^oo) with the 1 PN 
3-graviton vertex integral /oooooo- Also, there are 4 dis- 
tinct permutations, with the world line operator with 
2 graviton fields associated with either a, 6, c, or e. 
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= / dt 



r M a M b M c M e 



64(<i - 2)M 



pi 



- d ~ 1 

TT 2 



x (R-ae ^000000 c, e] + i? 3 e d /oooooo[a, e] 
+ R-be ^000000 [a, c, e]) 



dt 
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1024 7 ri( d - 1 )Mp 1 ( ^ 2) (d-2) 3 



f n3- d r>3 — d 
I -"he K be 



d3 — d d3 — d 
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{ K ab H ae 



3 — d p3— d 
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e>3 — d d3 — d\ d3 — d 



p3 — d p3 — d\ p3 — d 
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(a) 4 permutations, a, b, c or e has 2 graviton fields. 
©00 00© 



©oo oo© 

(b) No permutations necessary. 
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(c) 6 permutations. See text for discussion. 

FIG. 14: 2 PN four body diagrams with graviton vertices 
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Fig. (|14| c) Referring to the discussion under the 
2 body counterpart of Fig. (fl"4"] c). the 6 permutations 
of particle labels are: 2 x (ab\ce), 2 x (ac\be), and 2 x 
(ae\bc). 

The master integrals for Fig. (fj"4"] c) can be found 
in (|34l35l36l37p . -^oooo-oooo i s the linear combination 
of products of (|B3[) . As we did in the 2 body case, 
Iqqoo_oooo can he done using (|B3j) . its derivatives, and 
derivatives on 13. We have given the explicit expres- 
sion for /qooo- 0000 m the 3 body case ([40|l . but the 4 
body one is too lengthy to display. 

-^0000-0000 [*?' r ' s ' u ] i s -^ e ^ f° r possible future work. 

In terms of the ioooo-0000 m (|34I35I36I37|) . we have 
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Fig. (ll4| b^ Another straightforward application 
of Iqooooooo in (|33|) and (|B3[) provides us with 
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^. C[f 4 ] Effective Lagrangian 

Adding all the relevant diagrams, their permutations 
and the second order relativistic correction to kinetic 
energy from the rj^v^v" in the infinitesimal proper 
time ds now gives us the effective lagrangian describing 
the gravitational dynamics of n point masses at 0[v ] 
relative to Newtonian gravity in 3+1 dimensions: 



(2 PN) d=i ^2 Body j_ T 3 Body _,_ j 4 Body 



J cff 



(41) 
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2 Body 



2 2. \ 16 W - + 16 Wb 



G N M a M b 
Rab 



7 . 



Rab ■ V a -V a ■ V b - -V b ■ Vb } + Rba ' V b -V b ■ V a - -V a ■ V 



7 



I I Rab ■ Va 
8 I Rab 

3 

+ 7("!+«t- 2 i' v b ) 



Rba-V b +V b ) - - 



1 / Rba ■ V b 



Rab 



Rab ■ V a 



Rab ■ Va Rba ' Vb 3 (R ab ■ V a ) 2 (R ba ■ V b f 
Rab Rab 8 



Kb 



+ g (Rba ■ Vb vl+ Rab ■ V a vfj + ^R ab ■ V a R ba ■ V b 
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1 5 .!> 9 7 y , 1 , _» _ , \ 2 ^ -»2 — »2 ^ / — *2 -*2 \ — * — » 

+ yVa • % R a b + gK + "6 ) + ' «b) + ~V a V b - - (<U Q + V b )v a ■ V b 



G*M a M b ( 3 M a (Rba ■ v b ) 2 + M b {R ab ■ v a f 



B 2 

^ab 



B 2 

^ab 



2{M a + M b ) 



Rab ■ V a Rba ■ Vb 
^ab 



(2M a + M b )R ab ■ v a - (2M b + M a )B ba ■ v b 

v 2 a [2M a + Hm 6 ) + v% [2M b + Y Ma ) ~ ■ Vb (M a + M b ) 



GlM a M b ( 3, , ,A1 

RT^ [ MaMb + 2 {Ma + Mb) ) ) 

| \GlM a M b M c 

1 < r\ hi i - <L ft ^ 



L 



3 Body 



(42) 



l<a,b,c<.n 
,b,c distinct 



-v a + 8v b -v c ) + 



1 



RabRbc \2 



-v b + 8v a ■ v c + 



1 



RacRbc \2 



-v c + 8v a ■ v b 



Rba ■ V b Rca ■ V c , -R a b ■ V a R c b ■ V c , -R a c ■ "Ca-Rbc ' 4 
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RabRa 



RabRbc 



RacRbc 



vi vi 



Rab + Rac + Rbc \ Rbc Rac Rab 

''Rac ■ v a B ba ■ v b + Bac ■ v a B ca ■ v c + B ba ■ v b B ca ■ V c + 2(B ca ■ v c ) 2 
7 5 



2 Rab ^ac 

1 



RabRa 



Rba ■ V b + ~V a ■ V b + ~V b 



1 



(Rab + Rac + Rbc) 2 \R 
1 



-^2" ( 4i? a fc ' VcRba ' V b + 8B ab • V a R a b ■ V c - 2(B ab ■ v a ) z - 2(B, 



ab 



Vc)' 



RabRac 

l 



4-R a c • V b Rba ■ V b - 8B ac ■ VbRba ' V c + 12B ac ■ V a Rba ' ^c - 12-R QC • 'O'a-Rba ' V h 



1 



-Rab + -Rac + -Rbc V -R' 



ab 



—3- 8i? afc • tT Q -R afc • u c + 4i? ah • v c B ba ■ v b - 2(B ab ■ v c ) - 2(B ab ■ v a ) 



1 

Rab 

G 3 N M a M b M c 



2v 2 a — 4v a ■ v c — 2B ab ■ v a J j + 5 other permutations of {a, b, c} 
{M a + M c )B 2 ab 2M b B ab 3M a 



R 2 B 3 

ll ac ll bc 



R B 3 B 3 
n ac ih bc n ab 



if 3 

M a + ~M C I + 5 other permutations of {a, b, c} 



RabR 2 

+ ^2 ( Ma/ 22 [a,a,b,c} + M b I 22 [b, b,a,c}+ MJ 22 [ c , c , a, 6] ) - 2 + ^ + ^ 



ll Body 4 E G 3 N MaM b M c M e 



(43) 



1 



1 



1 



1 



l<a,b,c,e<n 
a,b,c,e distinct 

f I 22 [a, 6, c, e] _ 

[ 87T 2 \R a bRacRae RbaRbcRbe RcaRcbRce ReaRebRe 

1 / -Rbc 2i?^ e 2B a b\ 

Bab + -Rac + -Rbc \ RabRacRae RabR be R ae J 

+ 23 other permutations of {a, b, c, e} 



(44) 
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I 22 [q, r, s, u] ee lhn J d^y j d 3 ~ 2e z 5*6™ ((9fi?^ +2e )(ar n i?-; +2e )(i?- 1+2£ )(a;i?- 1+2£ )(9^-^ 2e ) 

after integration 

+ 5 permutations of {x , x r , x s , x u \ 1 R qz = \x q — z\ (45) 



The permutations in the definitions of Body and 
h\ Body means one would have to take the terms in 
the given square brackets [...], consider the resulting 
expressions obtained from permuting the particle la- 
bels as stated, and sum them all up at the end. For 
l22[q, t, s, u], if (qr\su) represents the term with x q and 
x r occurring in the y integration and the x s and x u in 
the z integration, then the 6 permutations in the defi- 
nition of I22 are: 2x (qr\su), 2x (qs\ru) and 2x (qu\rs). 

Relation to Ladm As a (partial) check of these 
results, we shall construct here a coordinate transfor- 

(2 PN) 

mation that would bring the 2 body portion of L y cS ' 
into the 2 body, acceleration-independent, lagrangian 
£adm in the literature; for example, eq. (178) of 
Blanchet @. (This construction can be found in 
Damour and Schafer SHI.) First, we note that defin- 
ing 



xi[t]=zi[t]+szi[t], 

where 8z a is assumed to be small relative to z a (an 
assumption to be justified shortly), would modify the 
form of the lagrangian L[{x a , v a , ■ • • }] up to first order 
in Sz in the following manner: 



L[{x a ,v a ,v a , ...}} 



SL 



= L[{z a , z a , z a , . . . }] + — 5z + total derivative 

oz 



SL \ - ^ ( d 

- z bz ^ E E 



dL 



\dt s d(d s z a /dt s 

Ka<n s=0 v v a > 



■ 6z a 



In particular, varying the Newtonian lagrangian 
Lq pn gives us 



E 



5L pn 

Sx a 



Sx a 



E 

l<ij"<3 l<a,b<n 



E j 



E 

i,j;a,b 



Rt 



xi + ^^Rl ) 



R 3 



j G N M a j , , 



G N M b j G N M a uJ ~ 

— H3 -"-ah — ~ -rti,„-^ 



/? 3 



ba 



(47) 



Because the term in the first line on the 
right hand side of (|4"T|) contains the "square" 
of the Newtonian equations of motion, namely 
(l/M a )(SL PN /5zi)(l/M b )(SL pn/Sz{), and because 
(1/M)(5L pn/Sz + SLi pn/Sz + ...)= 0, we see that 
this first term on the right hand side of (|4"T|) scales 
as [{1/M)6L PN /6z} 2 ~ [{1/M)8L X PN /^] 2 ~ [v 4 /r] 2 
and therefore can be discarded at the 2 PN order. 

Removal of accelerations The reason for lin- 
earizing the acceleration dependent terms in L 2 pn, 
keeping only the second and third lines on the right 
hand side of (|47|) is this. Denoting this linearized form 
of L 2 pn as L l 2 PN and referring to the —M a x a ■ 8x a 
piece in (I46|) , we see that the remaining acceleration de- 
pendent terms in L 2 PN can now be removed by defin- 
ing 



5z® ee 1 dLl PN 



M a dx a 



(48) 



Further transformations One can perform 
further coordinate transformations without re- 
introducing acceleration dependent terms. The key 
is to make the —M a x a ■ 8x a piece in (|4*fJ|) part of a 
total time derivative. Observe that, by having some 
arbitrary functional F depend only on positions and 
velocities, 12 we have the identity 



J2 M ° ( 



G N M b 
RL 



Rab ■ Sx a 



(46) 



/ y x a ~L, [{ x a,XaY\ 



Before proceeding with any coordinate transforma- 
tion, however, one needs to first re-write the terms 



quadratic in accelerations, Y^ij-ab^a^l^ ^ 



(2 PN) 
eff 



a.s 



12 We also exclude the possibility that F depend on time explic- 
itly, since our post-Newtonian lagrangian does not. 



2G 



[{x a ,X a }\ - 



dt 



Therefore, by putting 



62W = -L™ 

M a 3z a 

we can replace - J2 a M a z a -8zi l) = - Z)„ z a -dF/dz a 

with 12a Za ■ dF/dz a . 

At this point, let us note that the alterations to the 
form of the lagrangian due to 5z (|4"B")) occurs solely 
at the 2 PN order: 5z ~ rv 4 . Hence Sz is indeed 
small relative to z a , and it is only necessary to con- 
sider SLqpn/Sz and not SLipjq/dz, 6L2pn/5z, nor 
any corrections that are quadratic or higher polynomi- 
als of 6z . 

Altogether, the 2 body lagrangian after linearizing 
the accelerations and after the transformation x a = 
, less total derivative terms, now reads 

L[{x a ,v a ,v a , . . . }] 
= L pn[{z q , z a }} + Li PN [{4, z a }} 

+ L 2 p^\{z a , z ai z a}] 



E 
E 

a^b 



G N M h 

n-ab 



-Ral 



dL 



2 PN 



dx a 



dF 

dx a 



G^M b 

n-ab 



Rab 



8F \ 
dx a J 



(49) 



(50) 



It remains to construct F[{x a ,x a }]. From 
(l/M)(dF/dx) = 5z ~ rv 4 , we must have F ~ Mrv 5 . 
From dSDJ), we also must have F oc M a M b jM 2 x ~ 

G^M a M b ~ Mrv 2 , since all terms in L^ S PNS> need 
to contain at least one power of the mass of each of 
the 2 point particles and at least one power of Gn 
is required because all the terms (less the M-i^/16) 

(2 PN) 

in L eS are at least linear in Gn- To supply the 
additional v 3 needed, we have to consider all possi- 
ble products of the dimensionless scalars built out of 
terms occurring at 2 PN order, namely, {v 2 ,v 2 ,v a ■ 
v b , G N M a /R ab , G N M b /R ab } for the v 2 terms and {v a ■ 
R a b/Rab,v b ■ Rba/Rab} for the v 1 terms. The most 
general F is thus 



F[{£ s ,x s \s = l,2,...,n}] 
= 1 ^ / GlM a M b { Cl M a - 



a=£b 

G N M a M b 

Rab 



V a ■ Rab 



V a ■ Rab C 3 V a + C 4 V b + C 5 V a ■ V b 




+ (a <-> 6) 



(51) 



where the {ct|i = l,2,...,8} are arbitrary real num- 
bers, and (a <-> b) means one would have to take the 
terms occurring before it and swap all the particle la- 
bels a <-*■ b. 

Computing (|49|) and ([50]) with such a F reveals that 
one would recover the Ladm in Blanchet [|| from (|42[) 
for 



ci = 0, c 2 



3 n 1 

- — , C-i = U, Ca = — 

4 2 

cs =0, c 6 = 0, c 7 = -c 8 . 



where cs can be an arbitrary real number. 



IV. SUMMARY AND DISCUSSION 

In this paper we have, following (l3| . used a point 
mass approximation for the n-body in general relativ- 
ity, allowing us to obtain a lagrangian description at 
the cost of introducing an infinite number of terms 
in the action. Because we are seeking the 2 PN ef- 
fective lagrangian, however, only the minimal terms 
{— M a J ds a } are necessary. By examining the phys- 
ical scales in the problem, we have described how to 
organize our action ^ and outlined an algorithm that 
would allow us in principle to generate all the necessary 
Feynman diagrams up to an arbitrary PN order for a 
given set of point particle actions (minimal or not); 
as well as automate the computation process so that 
as few of the Feynman diagrams as possible are left 
for human evaluation. This way, the post-Newtonian 
program can be pursued in an efficient and systematic 
manner within the framework of perturbative held the- 
ory, and the necessary software may be developed to 
tackle the effective 2 body lagrangian calculation at 4 
PN and beyond. In the bulk of this work, we obtained 
in closed form the conservative portion of the effective 
lagrangian L c ft[{x a , v a , v a }] up to 1 PN for the general 
case of n point masses in d > 4 spacetime dimensions 
and up to 2 PN for 2 point masses in (3+l)-dimensions. 

It is apparent that the primary bottleneck of higher 
post-Newtonian calculations is one of calculus. For the 
n-body problem, it is the analytic evaluation of inte- 
grals such as I 22 (HSJ) and the 7i...jv (|B8|) for A > 4. 
At the same time, it is possible that choosing a dif- 
ferent gauge from the one used in (J7J) and/or a dif- 
ferent parametrization of h^ u may help reduce the 
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number of diagrams and the amount of work needed 
in manipulating the momentum dot products from 
the tensor contraction of graviton vertices in fourier 
space. The reduction of diagrams at 2 PN was re- 
cently demonstrated in the 2 body case computed by 
Gilmore and Ross 

0]. Th ey used the full de Don- 
dcr gauge 13 S gf = J d d Xy /gg^ 'g^ 'g<">T ' pa p and 
the Kaluza-Klein parametrization for first intro- 
duced to the PN problem by Kol and Smolkin [lj|. 14 
Some other possible choices include the ADM variables 
normally associated with the (3+l)-decomposition of 
the spacetime metric. Yet another possibility is to em- 
ploy the gravitational lagrangian constructed by Bern 
and Grant using quantum chromo dynamics gluon 
amplitudes, up to the 5-graviton interaction; this is 
sufficient, however, only to 3 PN. 

We end with a cautionary remark against taking 
superposition too literally within the post-Newtonian 
framework. By using the 1 PN lagrangian |23|) in 3 
spatial dimensions, and taking the continuum limit, 
the force F experienced by a stationary point mass M x 
located at r away from the center of a static, spherical, 
hollow shell of surface mass density a and coordinate 
radius R can be shown to be 15 



F = M. 



d r 



-Mi 



d 



1 + 12nRG N a dr \ \r\ 



In 



\r\+R 



\r\-R 



which evidently diverges as the point mass ap- 
proaches the surface of the shell. Because the force 
would vanish if gravity were purely Newtonian, such a 
result for a first order calculation most likely indicates 
the breakdown of perturbation theory in this regime, 
since the post-Newtonian lagrangian was derived with 
an implicit assumption that the point masses involved 
were well separated, i.e. r s <C r. 
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APPENDIX A: THE AT— GRAVITON FEYNMAN 
RULE 



Here we outline an algorithm that can be imple- 
mented on symbolic and tensor manipulation software 
such as Mathematica [2(| and the package FeynCalc 
[2~D | , to generate the Feynman rule for the N > 2 gravi- 
ton vertex in Minkowski space. 

Given a product of two function(al)s the 
term that contains exactly n powers of h is a discrete 
convolution 



(f-g\n)[h] = ^(f\m)[h](g\n-m)[h] 

m=0 

where (A|m)[/i] denotes the term in A that contains 
exactly m powers of h. Here we also assume that both 
/ and g can be developed as power series expansions 
starting from the zeroth power in h. 

With this observation, the term in the Einstein- 
Hilbert lagrangian containing exactly n powers of the 
graviton field is given by 



3 This will modify the 7V-graviton interaction in the Einstein- 
Hilbert action to all orders in h^ v . 

14 Comparing the calculus involved, however, when using the 
Kol-Smolkin parametrization for h^y, one notes that Gilmore 
and Ross for the 2 body problem, encountered the same 
master integrals as the ones used in this paper, dB3t and (IB4b ■ 
Moreover, for the n-body case, the most difficult integrals at 
2 PN, arising from the contraction of two 3-graviton vertices, 
namely ^Jooo-OOOO ll' r ' s < u l m GZJi arc identical in form to 
the ones in their equation (52), keeping all particles distinct. 

15 This computation arose out of discussions on Birkhoff's the- 
orem in GR with Dai De-Chang and Glenn Starkman. The 
following formula can also be found in their recent paper with 
Matsuo 0| • Its derivation actually only requires the 2 body 
portion of the L e g, because the 3-body portion integrates to 
a constant. 



-2M 



where 16 



(-),,, 



lb n 

' / d d x (Vlffl m) (R\n-m) (Al) 

m=0 J 



(A2) 



We are absorbing the ^Z 2 ' i n to the to save clutter. 
The full JV graviton rule would therefore be multiplied by a 



factor of M 



would contain no M„ 



JV(l-(d/2». 
pi 



for instance, the 2-graviton "vertex" 
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(^I n ) = ^ 



exp 



1 " 

IE 



s=l 

Tr( V - 1 hy = h v ' vl tr va ...h t '-\ 

n 

(n\n)=J2(9\mf v (R\n-mr 0av 



rn=0 



aB _ I '? 



(ffl») ap = 



n af3 



if n = 



(-)"^ q Mi ^m 2 ■■■h^- lP ifn>0 



Ml M 
,n.\ 



(i?|n) a ^ = 9 M (r|n)V-(M^^) 

n 

+ £ (r|m)% A (r|n - m) V - 0* <- 

m=0 



The complete Feynman rule for the N graviton ver- 
tex would be found by summing up the results from 
the above procedure for all the N\ possible permuta- 
tions of the external indices. An example featuring the 
3-graviton Feynman rule can be found in appendix B 

of m. 



APPENDIX B: INTEGRALS 

In this section we review the techniques involved in 
performing the Feynman integrals encountered in the 
n body problem at 2 PN. 18 

The starting point is the observation (usually at- 
tributed to Schwinger) that one may employ the inte- 
gral representation of the Gamma function, r[z]/& z = 
/o°° t z ~ le ~ M dt, for Re [b] > 0, to first give us the for- 
mula for combining multiple denominators. 



Note that the action for GR contains a total d- 
derivative term -2M^f 2 / d d x rf v (d a T% v - d^ a ), 
which needs to be discarded when deriving the Feyn- 
man rules. Furthermore, these Feynman rules will also 
be modified accordingly when a gauge fixing term is 
added. 

Because the graviton field is symmetric in its indices, 
to obtain the N graviton vertex with external indices 
{ai,/3i}, . . . , {on,Pn} given the action (|A1|) contain- 
ing exactly N powers of the graviton field, we first 
choose one particular set of contractions between the 
graviton fields in (jAljl with the N external ones. We 
then replace each field in (|A1[) with the identity 

tensor I" 4 ' 3 * carrying the indices {ai,fy} that cor- 
respond to those on the £th external field it is 

. ... a iPe 

being contracted with. The identity tensor reads 



da s a° 5 



1 

1 - Er=l 



<7\-\ hffJV 



(Bl) 



and second, together with the gaussian integral 
J exp[— x 2 ]dx = 7T 1 / 2 , to further yield 



and 



d A z _ n^T [a - |] 



(z2 + A) CT r[a]A CT -i 



(B2) 



I 



<*ePe 



2 (/V 



+ 



In momentum (k— )space, if we define the direction 
of momentum to be always flowing into the N gravi- 
ton vertex, we would also replace partial derivatives 
occurring in (jAlj) using the prescription 



[if ■ 



where Xe is the Ath component of the d-vector k of 
the £th external graviton h at p t that is contracted with 
h 17 



dp e l P' s 
(2tt) a [p>]° 



r 



4 CT 7r A /2r[cr]|f| 
0, 



A-2c 



for x ^ 

for x = 0, 

(B3) 



where the second equality is to be understood within 
the framework of dimensional rcgularization. 
An important corollary of (|B1[) and (|B2|) is 



d A ; 



The sign in front of the momentum vector, i.e. —ik\ ( . . . 
vs. +ik\ ( . . . , is actually immaterial because every term in 
the Einstein-Hilbert action contains two derivatives; what is 



important is to maintain a consistent sign convention for the 
arguments of the exponentials, either exp[ip°a; — ip ■ x\ or 
exp[— ip°x° + ip- x], in the fourier transforms. 
18 A comprehensive textbook on evaluating Feynman integrals 
is Smirnov 1131 . 
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n x/2 T[ x _ ai]T[ x _ a2]T[ai 



02 — 



r> 1 ]r[ ( 7 2 ]r[A-a 1 -a 2 ][2? 2 r+ 



(B4) 



By considering single and double spatial derivatives 
of (|B4[) . we may also obtain the formulas: 



T[pi H h pn]S 1 - X^Li a 



(z 2 + J2r. s =i (a r a s x^ - a r a s x r ■ x s ) 



piH hpw 



where the constraint ^ a r — 1 as well as a shift in 
the variable z have been employed. Writing 



(x - zf 



\z\Pi\x-z\Pz 



(B5) 



7T A /2 r 


A-pi 
2 


r 


A-p 2 - 


-2 


r 


P1+P2-A 
2 


a; 4 


r 




[f 


r[A 




P1+P2-2T 


|f|Pi+P2-A 



^A/2 r r [ r 1+ r 2 -2-A ] p [ 



(B6) 



A— r 2 +2 j 



2r ft] r [f ] r [a 

x ^(A-t 2 + 2)(ti+t 2 -2 
A-n 



n+T 2 -2 j 



(2(A + l)-(n+r 2 )) 



1ri+T 2 -A 



1T1+T2-2-A 



d A^ Oa - - 

\z-X a \P-\z-X b \P" 
7T A /2 r 



(B7) 



"A-p +2~ 


r 


"A-p b +2 


r 


"p a -2+p b -A 


2 




2 




2 



2r [fe]r[f]r[A- 



Pa+pb-4 j 



|Pa+P6-2-A 



(Pa + Ph - A - 2) 



\Xb - X a \ p " +Pb - X 



The last formula (|B7|) has been derived from (|B4[) by 
re-defining x = x a — x^ and shifting integration vari- 
ables before performing the appropriate derivatives. 

N point integrals We next review the evaluation 
of the "TV-point" integrals first carried out by Boos and 
Davydychev HEH. 



h.. 



N 



l — 



d x z 



U\(x s -z)^ 



(B8) 



These integrals can be viewed as the higher TV gen- 
eralizations of the TV = 2 case in (|B4[) . 
Applying (|B1[) transforms it into 



I\...n = 




N 

E 



1 N 

- 22 a rRrs 



r,s—l 
Rrs = l^r 



and recalling (|B2|) then allow us to deduce 



da s ag s 



[2 X/q ^p^g^pg]^*" ' 



(B9) 



By viewing a p R pq a q as components of a iV x iV sym- 
metric matrix with zeros on the diagonal, one notes 
that there are TV (TV — l)/2 distinct terms in the sum 
in the denominator of (|B9|) . To make further progress, 
one needs to iterate (L — 2) times the Mcllin-Barnes 
(MB) integral representation 



1 1 



(X + Y) T T[t] 2m J_ 



Y s 



dsr[ s + r]rH— , 



to obtain the [L — l)-fold MB representation for the 
sum of L terms in a denominator raised to some power, 



1 



1 



1 



\L-1 




[ / d Sa ri- 
x r [si + s 2 

In these integrals, the contour for the ith variable Sj 
is chosen such that the poles of the Gamma functions 
of the form r[- ■ ■ — Sj] lie to the right and those of the 
Gamma functions of the form r[- ■ • + Sj] lie to the left. 

Using (|BT0|) on (JB9J), followed by (JET]) with A a = 
= • • • = A-n(n—x)/2 = 1 and some careful algebraic 
reasoning, one arrives at the final form of the MB rep- 
resentation for the TV-point integral, 
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d z 



i 



(AT/2)(7V-1)-1 

nf^-W "r[A-E rfr ]nf =1 rWM M(f:i511 { ,JELa 

of N X N matrix 



dSy-r[-Sy] 



x r 



x r 



(iV/2)(iV-l)-l 



N 



{k.l}Guppcr A 
^ N 



n r 

l<s<iV 



Ps 



E 



S,'fc 



{j J fe}£5 t row and s* column 
upper A 



E + ft' - E 



{fc,Z}£uppcr A 



N 



J2 Pr + Pi' - 



Ski 



r=l 



{fc,/}(Eupper A 



(Bll) 



Here, {«',.?'} is some fixed pair of numbers chosen 
from the "upper triangular" portion of the N x N ma- 
trix of number pairs corresponding to their coordinates 
on the matrix; namely, the first row on the "upper A" 
reads from left to right, {1, 2}, {1, 3}, . . . , {1, N}, the 
second reads {2, 3}, {2, 4}, . . . , {2, N}, and so on un- 



til the (N — l)th row, which has only one element, 
{N-1,N}. 

N = 3 The MB integrals for the N = 3 case have 
been explicitly evaluated by Boos and Davydychev Q . 
One has from (jBlip . 



h 



23 



3 r tz>2 
l K 2 



r[A-£ r fr] n r r[ Pr ] (27 



du 



x r 



U + V + 



r [ Pl + u + v] r 



Pi - P3 



-^12" 








A 




2 ~ 


Pl - P2 



(B12) 



Assuming there exists a series representation of the 
integral /123 in powers of R 12 and -R13, we close the 
u- and v-contours on the right, turning each integral 
into an infinite sum over its residues by noting that 
r[— z] has singularities on the complex z plane only 
in the form of simple poles at zero and the positive 
integers. (Whether one should close the contour to 
the left or to the right really depends on the numerical 
range of R12 and R13 considered. For instance, the MB 
representation (1 + z)~ x = (2mT[\\)- 1 f+*™ duT[u + 
A]r[— u]z u can be converted into a power series in 1/z 
or z, for A > 0, by closing the contour to the left or 
right, depending on whether \z\ > 1 or \z\ < 1. Here 
we will simply assume, for each choice (left or right), 
there is some region of i?i2,i?i3 £ 1 in which it is 
valid.) The residues of r[— z] at these locations are 



Res T[-z] 



m = 0,+l,+2, . 



r[- u] and 2 of the form r[ v], this converts the 

2-fold MB integrals into a 2-fold infinite sum of 4 terms. 
One can proceed to change summation variables and 
manipulate the Gamma functions in the summands us- 
ing the definitions for the Gauss hypcrgeometric func- 
tion and the Pochhammer symbol 



2 F 1 [a,b;c;z] = ^ 



u\(c) u 



(a) u = a(a + 1) . . . (a + (u — 1)) 



T[a + u] 



and the relation 



Because there are 2 Gamma functions of the form 



(-)"TH 

r [t - m = V li, meZ,reC 

(1 - T) m 

to further reduce the 2-fold sum into a single sum: 
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2^1 
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2 2 Kb 



(B13) 



As described in Boos and Davydychev this sum 
has a closed form expression in terms of the Appell 
hypergeometric function F4 of two variables, which has 

a perturbative definition of the form xJ ( a )jWj 77 r * 1 • x 1 

2^7i M . 2*i[-J.i-7-j;*;y] 

(a) m+n (/3) m+n y" = ^4 [a, 13; 7, «5; s, zy] , 



F 4 [a, /3; 7, 5; x, y] = 2J X! 



Through the relation we now have 



Iabc[PaiPb, Pc 
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[(J?- £„) 2 ]'»[(z - ?i) a ]«[(^- ^c) 2 ]^ 
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However, F4 is not defined in Mathematica [2(j |. 
whereas the sum in (|B13|) can easily be entered. In- 
particular, at 2 PN, the n body problem requires the 
knowledge of 



6 ij d?d?I 123 \p a =p b = Pc = 1/2] = -4ir(R ab R a 



h[a,b,c] 



4:77^ 



'123 



rf-3 



Pi = P2 = P3 



Applying, in 3 — 2e spatial dimensions, the Laurent 
expansion for the Gamma function about negative in- 
tegers or zero, 



(-) m 1 



m 1 \ 

r=l / 



m = 0,1,2,..., 

to the summand in (|B13[) . before employing the 
Mathematica command FullSimplify on the summation 
(IB13|) . yields the final form for I 3 [a, b, c], 



I 3 [a,b,c] 
1 



64tt 2 



- + 2 - 27e - 2 ln[7r] 



-4ln\R ab + R ac + R bc \+0[e]j (B14) 

where 7e = 0.57721 ... is the Euler-Mascheroni con- 
stant and the hyperbolic function identity tanh _1 [z] = 
(1/2) (In |1 + z\ — hi 1 1 — z\) was used. An alternate 
derivation of this result can be found in Blanchct et 

al. [n|. 

A direct computation would show that this result 
is consistent with the Poisson equation obeyed by the 
A = 3 integral in 3 spatial dimensions, 



APPENDIX C: 3 PN DIAGRAMS 

In this section, we collect the fully distinct Feynman 
diagrams necessary for the computation of the effec- 
tive lagrangian for n non-rotating, structure-less point 
masses as described by the minimal action in ^ at the 

3 PN order. Fully distinct here means that, to obtain 
the full 3 PN lagrangian one would have to, whenever 
applicable: 

• Consider all possible permutations of the particle 
labels of the diagrams displayed. 

• For the n = 3,4 and 5 diagrams, consider all 
possible ways of setting some of the particle la- 
bels equal to each other, so that from the n = 3 
diagrams one would obtain their n = 2 counter- 
parts; from the n = 4 their n = 2 and 3 counter- 
parts; and from the n = 5 their n = 2,3 and 4 
counterparts. 

The 2 body diagrams are in Fig. (fT5|) . The 3 body 
diagrams with graviton vertices are Fig. (fT!)l [T7|) ; and 
those with no graviton vertices are Fig. (|T5]) . The 

4 body diagrams with graviton vertices are Fig. (|19[ 
l20l l2"Tj) ; and those without graviton vertices are Fig. 
(|2"2l |2"3")) . Finally the 5 body diagrams with graviton 
vertices can be found in Fig. (j2"5|) , whereas those with 
none can be found in Fig. (|24|) . 



33 



2 

^00 00 ^ ^00 v.OO^ ^Oi 0/>^ ^ij 00^ 

0--X--X--X-® 0- -x X-® ®---x x - - ® --x x--® 

2 2 4 2 2 

^ 00, v v , 00^ ^ 00 v, 00^ ^ Oi v, O/V-x a-n 00 v , ij ^ 

0------X------® x ® © X-------0 x----- y g 

2 2 4 2 2 

®ij v , 00 1 — . / — N // //' 00, v v , 00 s~ N 0/, v \/ , 0/ 
X--- -® ® y x y © • - •■■© ©-- ® 

2 2 6 4 4 

v v , 00 00 v , 00^^ 00 v , ij ij v , 00^ 

© © © © ® 
4 2 

FIG. 15: 3 PN 2 body diagrams. 
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FIG. 16: 3 PN 3 body diagrams containing 3 graviton vertices, 1 of 2. 
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FIG. 17: 3 PN 3 body diagrams containing 3 graviton vertices, 2 of 2. 
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FIG. 18: 3 PN 3 body diagrams with no graviton vertices. 
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FIG. 19: 3 PN 4 body diagr 
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with graviton vertices, 1 of 3. 




FIG. 20: 3 PN 4 body diagrams with graviton vertices, 2 of 3. 




FIG. 21: 3 PN 4 body diagrams with graviton vertices, 3 of 3. 




FIG. 24: 3 PN 5 body diagrams with no graviton vertices. 




FIG. 25: 3 PN 5 body diagrams with graviton vertices. 



